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Multiply _ 

feet 

inches 

gallons  (US  liquid)  per  minute 


-M _ 

0.3048 

2.54 

0.000006309 


SI  (metric)  units  as 

_ To  Obtain _ 
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LEVEEMSU :  A  SOFTWARE  PACKAGE  DESIGNED  FOR 
LEVEE  UNDERSEEPAGE  ANALYSIS 


PART  I :  INTRODUCTION 

Background 

1.  A  Repair,  Evaluation,  Maintenance,  and  Rehabilitation  (REMR)  Levee 
Underseepage  Workshop  was  held  at  the  US  Army  Engineer  Waterways  Experiment 
Station  (WES)  on  10  April  1984  to  establish  research  needs  related  to  levee 
underseepage  control.  Representatives  from  the  Rock  Island,  St.  Louis,  Mem¬ 
phis,  and  Vicksburg  Corps  of  Engineers  Districts  attended.  One  research  task 
identified  was  comparing  predicted  levee  underseepage  conditions  to  observed 
performance.  This  task  emerged  because  of  concerns  that  Corps'  procedures  and 
criteria  may  be  overly  conservative  in  many  instances,  requiring  costly  con¬ 
trol  measures  where  they  may  not  be  needed,  but  may  be  unconservative  in  other 
cases  by  failing  to  identify  areas  where  sand  boils  may  occur. 

2.  In  September  1986,  a  critical  review  of  underseepage  analysis  proce¬ 
dures  was  prepared  by  Wolff  (1986).  This  review  noted  that  Corps'  analysis 
and  design  procedures  required  a  high  level  of  judgment  in  formulating  the 
problem  for  analysis.  In  particular,  actual  soil  profiles  and  topography  are 
often  quite  irregular,  but  the  available  procedures  required  modeling  level 
topography  with  uniformly  thick  soil  layers.  The  judgments  required  for  this 
step  alone  could  result  in  greatly  different  analyses  by  different  designers. 
The  study  reported  herein  addressed  the  development,  testing,  and  use  of  a  new 
computer  program  LEVEEMSU  for  analysis  of  levee  underseepage  for  cross  sec¬ 
tions  with  irregular  geometry. 

Previous  Studies 


3.  The  Corps'  procedures  are  based  on  closed- form  solutions  for  differ¬ 
ential  equations  of  seepage  flow  presented  by  Bennett  (1946) .  For  conditions 
of  irregular  geometry  of  variable  properties,  sc1utions  cannot  be  obtained  in 
closed-form,  but  can  be  obtained  numerically.  In  1987,  initial  research  was 
conducted  at  Michigan  State  University  regarding  the  application  of  numerical 
methods  to  levee  underseepage  analysis  (Wolff  1987) .  It  was  shown  that 
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special-purpose  computer  programs  had  certain  advantages  over  both  traditional 
undeiseepage  analysis  procedures  and  general-purpose  seepage  analysis  pro¬ 
grams.  As  previously  noted,  traditional  procedures  (US  Army  Engineer  Water¬ 
ways  Experiment  Station  1956a,  1956b)  require  that  single  values  be  assigned 
to  the  variables  even  though  the  stratum  thicknesses,  ground  and  water  eleva¬ 
tions,  etc.,  often  assume  different  values  at  different  points  in  the  cross 
section.  General-purpose  seepage  analysis  programs  using  the  finite  element 
method  (e.g.,  Tracy  1973)  can  model  such  irregularities;  however,  they  often 
require  a  relatively  high  degree  of  effort  to  model  a  problem  and  interpret 
the  results,  even  when  pre-processors  and  post-processors  are  used. 

4.  The  1987  research  included  the  development  of  three  FORTRAN  codes: 
LEVEEIRR,  to  model  irregular  geometry;  LEVEE3L,  to  model  three-layer  founda¬ 
tions;  and  LEVEECOR,  to  model  corners  or  bends  in  levee  alignment.  These  were 
"preliminary"  programs  developed  to  demonstrate  the  feasibility  of  the  numeri¬ 
cal  approach.  These  programs  were  used  to  analyze  actual  data  at  a  number  of 
levee  reaches  and  back-calculate  field  permeability  values. 

Scope 

5.  This  report  documents  the  development,  testing,  and  use  of  a  new 
computer  program,  LEVEEMSU,  for  analysis  of  levee  underseepage.  This  program 
represents  a  second-generation  version  of  LEVEEIRR  described  above,  and  in¬ 
cludes  a  number  of  enhancements.  The  program  uses  numerical  methods  to  ana¬ 
lyze  underseepage  for  two-dimensional  levee  cross-sections  having  nonuniform 
geometry  and  properties.  Thicknesses  and  elevations  of  soil  layers,  the 
ground  surface,  and  ponded  water  all  may  vary  in  the  horizontal  direction. 

Top  blanket  permeabilities  may  be  specified  independently  for  each  side  of  the 
levee  and  may  be  constant  or  may  vary  as  a  function  of  blanket  thickness. 

Heads  and  gradients  are  calculated  as  a  function  of  horizontal  location.  The 

effects  of  a  line  of  relief  wells  may  be  modeled.  The  program  features  a 
graphic  display  of  input  and  results  to  aid  in  checking  the  input  and  inter¬ 
preting  the  results.  The  graphic  window  may  be  changed  by  the  user  to  look  at 
various  regions  of  the  solution  at  any  desired  scale.  The  program  is  particu¬ 
larly  useful  for  analyzing  and  designing  ditches,  borrow  pits,  etc. 
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6.  LEVEEMSU  provides  the  user  a  number  of  advantages  over  other  methods 
of  analysis.  As  the  analysis  of  interest  always  involves  a  levee  and  two  soil 
layers,  data  entry  can  be  made  more  concise  than  for  a  general-purpose  finite 
element  program,  and  the  program  can  be  designed  compactly  to  provide  rapid 
solutions  on  inexpensive  hardware  with  minimal  memory.  Likewise,  output  can 
be  arranged  to  provide  results  in  the  most  meaningful  form  (e.g.,  gradient 
through  the  top  blanket  versus  distance). 
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PART  II:  PROGRAM  DESCRIPTION  AND  SOLUTION  TECHNIQUES 
Program  Description 

7 .  The  computer  program  LEVEEMSU  is  an  entirely  new  program  based  on  a 
previous  program  named  LEVEEIRR  (Wolff  1987)  .  LEVEEMSU  was  developed  for 
analysis  of  levee  underseepage  and  design  of  underseepage  control  measures 
where  it  is  desired  to  model  cross  sections  having  nonuniform  geometry  or 
properties . 

8.  The  program  is  furnished  as  a  binary  executable  file,  LEVEEMSU.EXE, 
designed  to  run  on  IBM  (TM)  and  compatible  personal  computers  under  the  MS  DOS 
operating  system.  A  math  coprocessor  is  highly  recommended.  No  computer 
language  or  compiler  need  be  installed  on  the  computer.  The  program  was  de¬ 
veloped  using  Microsoft  QuikBasic  (TM)  and  linked  to  required  library  files  to 
produce  a  single  executable  file.  The  QuikBasic  language  was  selected  in  lieu 
of  the  more  traditional  FORTRAN  to  maximize  the  use  of  color  and  graphics 
capabilities  of  microcomputers  yet  retain  the  mathematics  of  the  source  code 
in  a  form  that  is  reasonably  readable  to  engineer  programmers.  The  program 
can  be  run  in  three  graphics  modes,  EGA  color,  EGA  monochrome,  and  CGA  mono¬ 
chrome  ,  depending  on  the  available  graphics  card,  monitor,  and  whether  a 
graphics  screen  copy  is  desired.  In  the  EGA  color  mode,  the  geometry  of  the 
substratum,  top  stratum,  water,  and  piezometric  grade  line  are  displayed  in 
color.  In  the  EGA  and  CGA  monochrome  mode,  these  are  displayed  in  high- 
resolution  and  medium- resolution  monochrome,  respectively.  In  the  monochrome 
modes,  the  graphic  screen  can  be  copied  to  a  graphics  printer  using  a  screen 
dump  program  such  as  GRAPHICS.COM  for  CGA  and  EPSON.COM  for  EGA. 

9.  The  program  reads  input  data  from  a  separate  data  file.  The  format 
of  the  input  file  is  described  in  Appendix  A.  The  program  displays  default 
values  for  certain  variables  which  affect  the  time  required  for  solution  and 
solution  accuracy.  These  values  can  be  changed  from  the  keyboard  during  pro¬ 
gram  execution.  Results  of  the  analysis  are  displayed  on  the  graphic  screen; 
a  detailed  summary  of  the  results  is  written  to  an  output  file  which  can  be 
printed  during  program  execution  or  separately  later.  Details  on  running  the 
program  are  described  in  Appendix  B.  An  example  run  is  shown  in  Appendix  C. 
Standard  input  data  files  are  discussed  in  Part  III  and  listed  in  Appendix  D. 
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Appendix  E  presents  a  hand  calculation.  A  program  listing  is  shown  in  Appen¬ 
dix  F. 


10.  Subsurface  conditions  beneath  levees  in  alluvial  valleys  are  tradi¬ 
tionally  modeled  as  two  soil  layers,  a  semipervious  top  blanket  or  top  stratum 
of  clay,  silt,  or  silty  sand  overlying  a  pervious  substratum  of  sand.  Flood 
conditions  riverside  of  the  levee  result  in  downward  flow  of  seepage  through 
the  riverside  top  blanket,  lateral  flow  through  the  pervious  substratum,  and 
upward  flow  through  the  landside  top  blanket.  Given  certain  conditions  of 
geometry  and  soil  properties,  the  upward  gradient  in  the  landside  top  blanket 
can  be  excessive,  and  safety  against  boiling  is  of  concern.  Underseepage 
analyses  are  performed  to  predict  the  piezometric  head  along  the  base  of  the 
landside  top  blanket  (or  at  least  at  the  levee  toe)  and  the  gradient  through 
the  blanket  as  functions  of  riverside  and  landside  water  levels.  Where  cal¬ 
culations  indicate,  excessive  gradients  are  expected,  and  control  measures  are 
designed.  These  are  typically  seepage  berms  or  relief  wells.  Additionally 
analyses  may  be  performed  to  assess  the  effect  of  proposed  or  existing  control 
measures . 

11.  A  solution  for  the  piezometric  head  beneath  a  semipervious  top  blan¬ 
ket  adjacent  to  a  dam  or  levee  on  a  pervious  substratum  was  proposed  by 
Bennett  (1946).  Bennett  assumed  perfectly  horizontal  flow  in  the  pervious 
substratum  and  perfectly  vertical  flow  in  the  top  blanket.  If  the  thicknesses 
and  permeabilities  of  the  blanket  and  the  substratum  are  taken  as  constants, 
the  piezometric  head  at  the  base  of  the  blanket  and  the  upward  gradient 
through  the  blanket  can  be  directly  calculated  for  a  number  of  various  boun¬ 
dary  conditions  using  equations.  Solutions  have  been  widely  published  within 
the  Corps  of  Engineers  (US  Army  Engineer  Waterways  Experiment  Station  1956a, 
1956b;  Office,  Chief  of  Engineers  1986a,  1986b)  and  elsewhere  (Turnbull  and 
Mansur  1961)  .  Underseepage  analysis  by  the  Corps  traditionally  has  utilized 
these  closed-form  solutions. 
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Analysis  of  Irregular  Geometry 

12.  If  numerical  methods  are  used  to  solve  Bennett's  (1946)  differential 
equation,  the  foundation  geometry  and  properties  need  not  be  uniform.  Rather, 
values  can  be  assigned  or  interpolated  at  a  number  of  points  or  nodes  (as  many 
as  desired) ,  and  the  differential  equation  satisfied  approximately  at  each 
node.  Solution  techniques  have  been  presented  by  Wolff  (1987)  and  are  extend¬ 
ed  herein. 

13.  A  unit  width  of  levee  is  modeled  as  a  two-dimensional  cross  section. 
Seepage  flow  is  assumed  to  be  horizontal  in  the  substratum  and  vertical  in  the 
top  blanket.  A  one -dimensional  numerical  solution  is  obtained  by  considering 
a  line  of  nodes  along  the  interface  between  the  substratum  and  blanket.  The 
program  user  describes  the  foundation  geometry  using  x  and  y  coordinates 
along  a  number  of  vertical  sections,  in  a  fashion  similar  to  the  data  input 
for  slope  stability  analysis  programs.  The  program  generates  a  set  of  nodes 
and  associated  geometry  information  based  on  the  user  input.  Dimensions  and 
properties  are  assumed  to  vary  linearly  between  nodes.  As  the  piezometric 
head  in  the  substratum  is  implied  to  be  constant  along  any  vertical  section, 
the  node  actually  represents  the  entire  thickness  of  the  substratum. 

14.  Figure  1  illustrates  conditions  at  a  typical  node.  The  node  (J)  is 
located  at  coordinates  XX(J)  and  YY2(J)  .  In  the  x  direction,  the  node 
represents  a  length  of  substratum  and  blanket  extending  halfway  to  each  ad¬ 
jacent  node,  XX(J-l)  and  XX(J+1)  .  In  the  y  direction,  the  node  is  asso¬ 
ciated  with  a  substratum  thickness  D(J)  -  YY2(J)  -  YY1(J)  ,  a  blanket  thick¬ 
ness  Z(J)  -  YY3(J)  -  YY2(J)  ,  and  a  landside  water  elevation  YYWATER(J)  . 

The  piezometric  elevation  at  the  node,  PIEZEL(J)  ,  is  calculated  by  the 
program. 

15.  Flow  through  the  represented  element  of  the  foundation  is  lumped  at 
the  node  for  analysis.  Referring  to  Figure  1,  continuity  requires  that  at 
each  node  landside  of  the  levee, 

Qin  -  Qout  +  Qup  (eq  1) 

where  Qin  is  the  flow  in  the  substratum  toward  the  node,  Qout  is  the  flow  in 
the  substratum  beyond  the  node ,  and  Qup  is  the  flow  or  seepage  through  the 
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PIEZEL(J) 

ir- . 

YYWATER(J) 


XX(J-I)  XX(J) 


Figure  1.  Geometry  at  a  typical  node 


top  blanket  in  the  vicinity  of  the  node.  On  the  riverside  of  the  levee,  the 
same  equation  is  used,  but  Qup  assumes  a  negative  value.  Between  the  land- 
side  and  riverside  levee  toes,  Qup  is  taken  as  zero.  From  Darcy's  law 

Q  -  kiA  (eq  2) 

where  Q  is  the  flow,  i  is  the  hydraulic  gradient,  and  A  is  the  area  nor¬ 
mal  to  the  flow.  The  flow  terms  can  be  approximated  numerically  as  follows: 


| PIEZEL(J -1) -PIEZEL(J) |  |D(J)  +  D(J-l) [ 
^  ;  |  XX(J)  -  XX(J-l)  ||  2  | 


(eq  3) 


„  | PIEZEL(J) -PIEZEL(J+1) |  | D( J+l)  +  D( J) j 

y°ut  ^  '  |  XX(J+1)  -  XX(J)  II  2  I 


(eq  4) 


<up 


(KB(J) ) 


|  PIEZEL(J) -YYWATER(J) |  |XX(J+1)  +XX(J-1)| 
I  Z(J)  ||  2  | 


(eq  5) 


where 

KF  is  the  horizontal  permeability  of  the  pervious  substratum 
KB(J)  is  the  vertical  permeability  of  the  blanket  at  node  J 
D(J)  is  the  thickness  of  the  pervious  substratum  at  node  J 
Z(J)  is  the  thickness  of  the  top  blanket  at  node  J 
XX(J)  is  the  horizontal  location  of  node  J 

YYWATER(J)  is  the  elevation  of  ponded  water  (or  ground  surface)  at 
node  J 

PIEZEL(J)  is  the  elevation  of  the  piezometric  surface  at  node  J 


Substituting  the  flow  Equations  (2  through  5)  into  the  continuity  Equation  1, 
the  piezometric  elevation  at  any  node  J  ,  PIEZEL(J)  ,  can  be  expressed  as 


PIEZEL(J) 


PIEZEL(J - l)*Cl (J)+PIEZEL(J+1)*C2 (J)+YYWATER( J)*C3 ( J) 
C1(J)  +  C2(J)  +  C3(J) 


(eq  6) 
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where 


(KF)*(D(J)+D(J-1) ) 
U  '  “  (XX(J)  -  XX(J - 1) )*2 


(eq  7) 


^o/ts  (KF)*(D(J+1)+D(J) ) 

'  ~  (XX(J+1)  -  XX(J)  )*2 


(eq  8) 


C3(J)  - 


(KB)*(XX(J+1) -XX(J-l) 
Z(J)*2 


(eq  9) 


To  obtain  a  solution,  Equation  6  is  solved  by  iteration. 


Variable  Node  Spacing 


16.  Nodes  are  generated  at  the  x  coordinates  specified  by  the  user  and 
at  a  number  of  intermediate  locations.  The  locations  of  all  generated  nodes 
are  shown  on  the  graphic  screen  and  listed  in  the  output  file.  The  number  of 
nodes  used  for  analysis  affects  both  solution  accuracy  and  solution  time.  To 
optimize  both  of  these  factors,  the  node  generating  algorithm  in  LEVEEMSU  pro¬ 
duces  nodes  at  a  variable  spacing.  Near  the  riverside  and  landside  levee 
toes,  where  gradients  are  the  highest  and  change  most  rapidly,  nodes  are  gene¬ 
rated  at  a  maximum  distance  of  25  ft  apart.  The  distance  between  nodes  is  a 
default  value  and  can  be  changed  by  the  user  during  program  execution.  At 
progressively  further  distances  landside  and  riverside  from  the  levee  toes, 
nodes  are  spaced  increasingly  further  apart.  This  technique  and  the  spacing 
ratios  set  in  the  program  have  been  found  to  provide  reasonable,  fast,  and 
consistent  solutions  with  relatively  few  nodes.  The  algorithm  used  produces 
much  more  consistent  results  than  the  scheme  used  in  LEVEEIRR,  the  predecessor 
program,  which  generated  a  fixed  number  of  nodes  between  each  specified 
section. 


Lands ide  Water  Elevation 


17.  The  landside  water  elevation  can  be  specified  independently  at  any 
location.  Flow  is  assumed  to  occur  vertically  through  the  top  blanket  driven 
by  a  head  equal  to  the  difference  of  the  piezometric  elevation  at  the  base  of 
the  blanket  and  the  specified  landside  water  elevation.  Where  landside  water 
elevation  is  above  the  ground,  consistent  values  should  be  specified  to  model 
the  water  surface.  Where  the  landside  ground  is  irregular,  specifying  the 
water  surface  coincident  with  the  ground  surface  will  model  water  rising  to 
the  surface  and  running  off.  Where  landside  swales  are  separated  by  relative¬ 
ly  high  ridges,  the  user  may  wish  to  specify  landside  water  surface  elevations 
lower  than  the  ground  surface  under  the  crowns  of  the  ridges.  At  the  user's 
option,  water  levels  in  swales  may  vary  from  swale  to  swale. 

Variable  Blanket  Permeability 

18.  The  in  situ  vertical  permeability  of  a  uniformly  thick  top  blanket 
during  flood  may  be  significantly  different  on  the  riverside  and  landside  of  a 
levee.  On  the  riverside,  downward  flow  may  enhance  siltation,  plugging  of 
cracks  and  defects,  etc.,  reducing  the  effective  permeability.  On  the  land- 
side,  upward  flow  may  tend  to  open  defects  in  the  blanket,  increasing  the 
permeability.  These  differences  may  be  modeled  with  LEVEEMSU  by  specifying 
different  permeability  values  for  the  riverside  and  landside.  When  solving 
Equation  6,  the  program  will  check  to  see  whether  a  node  is  riverside  or  land- 
side  of  the  levee  and  assign  the  appropriate  value  for  KB  . 

19.  A  further  refinement  allows  the  permeability  to  vary  inversely  as  a 
function  of  blanket  thickness.  For  levee  design  along  the  lower  Mississippi 
River,  permeability  values  are  often  assigned  using  curves  of  permeabil  ity 
versus  blanket  thickness.  This  practice  reflects  the  greater  probability  of 
blanket  defects  in  thin  blankets  versus  thick  blankets.  Figure  2  shows  the 
relationship  between  top  blanket  thickness  and  permeability  used  by  the  Lower 
Mississippi  Valley  Division  (LMVD) . 

20.  It  is  logical  to  extend  this  practice  to  the  case  of  blankets  of 
variable  thickness.  For  this  case  of  a  ditch  or  borrow  pit  cut  partly  through 
a  clay  top  blanket,  it  is  reasonable  to  expect  a  higher  permeability  in  the 
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Permeability,  k,  ft/min. 


Figure  2.  Blanket  permeability  versus  top  stratum  thickness 
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ditch  bottom  or  pit  bottom  than  in  the  thicker  adjacent  blanket.  Rather  than 
program  the  curves  from  Figure  2,  it  was  considered  more  practical  to  provide 
the  analyst  an  infinite  number  of  curves.  This  is  done  by  assuming  permeabil¬ 
ity  versus  blanket  thickness  functions  that  are  a  family  of  straight  lines  on 
semilog  paper  approximately  parallel  to  the  LMVD  curves.  These  functions  are 
generated  by  the  following  equation: 


KB  - 


K10 

exp( -0 . 065924* ( 10 -Z) ) 


(eq  10) 


where  KB  is  the  permeability  for  a  blanket  thickness  Z  and  K10  is  the 
permeability  of  a  10-ft- thick  blanket  of  a  given  material.  The  resulting 
functions  are  superimposed  on  Figure  2  for  selected  values  of  K10  .  When 
running  under  this  option,  the  user  specified  a  "curve  number"  or  a  value  for 
K10  and  the  program  will  calculate  the  corresponding  blanket  permeability 
value  at  each  node . 


Modeling  a  Line  of  Relief  Wells 

21.  Rigorous  analysis  and  design  of  relief  wells  is  a  three-dimensional, 
nonlinear  problem  which  is  beyond  the  scope  of  this  report  and  computer  pro¬ 
gram.  However,  LEVEEMSU  is  capable  of  approximately  assessing  the  effect  of  a 
line  of  relief  wells  using  an  option  to  specify  the  piezometric  elevation  at 
one  x  coordinate.  A  specified  piezometric  elevation  can  represent  the  aver¬ 
age  head  in  a  line  of  wells  (Hav)  .  When  this  option  is  used,  the  variable 
PIEZEL(J)  is  forced  to  the  assigned  value  at  the  node  closest  to  the  speci¬ 
fied  location  and  that  node  is  skipped  in  the  iterative  solution  process.  The 
program  will  then  calculate  the  flow  to  the  well  line  required  to  achieve  the 
specified  piezometric  elevation  as  follows: 

Qwell  ”  Qin  "  Qout  "  Qup  (eq  11) 

This  is  illustrated  in  Figure  3.  The  analyst  can  then  use  conventional  meth¬ 
ods  to  design  a  well  system  that  will  pass  a  flow  of  Qwell  under  the  speci¬ 
fied  head  conditions. 
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PIEZEL(J)  =  YWELL 
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XX(J)  •  XWELL 


Figure  3.  Analysis  at  a  relief  well  node 
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PART  III:  PROGRAM  TESTING  AND  PARAMETRIC  STUDIES 


Closure  Tolerances,  Node  Spacing, 
_ and  Number  of  Iterations 


22.  The  program  solves  different  equations  at  a  finite  number  of  points 
to  approximate  the  solution  of  a  differential  equation  over  a  continuous  do¬ 
main.  No  exact  solution  is  obtained;  rather,  continued  iterations  produce 
successively  more  accurate  solutions.  The  iteration  procedure  is  stopped  when 
the  maximum  residual  (change  in  calculated  piezometric  elevation  at  any  node 
between  successive  iterations)  is  less  than  a  specified  closure  tolerance. 

The  program  incorporates  default  values  for  the  closure  tolerance,  node  spac¬ 
ing  near  the  levee,  and  maximum  number  of  iterations.  These  can  be  changed  at 
the  user's  option  during  program  execution.  The  closure  tolerance  must  be 
significantly  less  than  the  desired  accuracy  of  the  solution,  as  small  changes 
from  one  iteration  to  the  next  may  accumulate  toward  the  "exact"  solution. 
Smaller  node  spacings  and  closure  tolerances  yield  more  accurate  solutions  but 
require  longer  times  to  run. 

23.  A  standard  input  file  named  DATACHK  (listed  in  Appendix  D)  was  used 
to  initially  check  program  solutions  and  evaluate  relationships  among  toler¬ 
ance,  node  spacing,  and  number  of  iterations.  This  input  file  models  a  levee 
with  a  uniform  substratum  and  top  blanket,  D  -  80  ft  and  Z  =  10  ft  ;  founda¬ 
tion  lengths  L]*  -  1,500  ft  ,  Lz  -  300  ft  ,  L3  -  3,200  ft  ;  and  a  permeabil¬ 
ity  ratio  kf/kb  -  1,000  .  The  final  screen  output  for  a  run  using  DATACHK 
and  default  options  is  shown  in  Figure  4.  As  show;.,  the  program  calculates  a 
maximum  residual  head  of  8.87  ft  and  a  maximum  gradient  of  0.89,  both  occurr¬ 
ing  at  the  landside  toe  (x  -  1,800  ft). 

24.  The  sensitivity  of  the  results  to  node  spacing  and  tolerance  was 
examined  by  systematically  adjusting  these  values  during  program  execution  and 
noting  changes  in  calculated  maximum  gradient  and  number  of  iterations  re¬ 
quired  for  solution.  Results  are  shown  in  Figures  5  and  6.  Also  shown  in 
Figure  5  is  the  theoretical  maximum  gradient  of  0.88  obtained  from  a  hand 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  G) . 
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Figure  5.  Maximum  gradient  versus  tolerance,  file  DATACHK 


10  ft,  110  nodes 
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Figure  6.  Number  of  iterations  versus  tolerance,  file  DATACHK 


check  (Appendix  E) .  In  Figure  5,  it  is  seen  that  a  node  spacing  of  25  ft  or 
less  and  a  tolerance  of  0.0005  ft  are  sufficient  to  calculate  a  gradient  with¬ 
in  a  few  hundredths  of  the  theoretical  value.  A  node  spacing  of  50  ft  was  too 
coarse  to  accurately  calculate  the  gradient  due  to  the  averaging  of  upward 
flow  conditions  near  the  toe  over  too  large  a  horizontal  distance.  Figure  6 
illustrates  the  number  of  iterations  required  to  converge  within  the  specified 
tolerance.  Decreasing  the  minimum  node  spacing  to  10  ft  increases  the  number 
of  nodes  from  46  to  110,  and  increases  the  number  of  iterations  from  a  few 
hundred  to  a  few  thousand.  This  produces  a  negligible  increase  in  accuracy 
(Figure  5)  but  a  relatively  large  increase  in  solution  time. 

25.  The  analyses  just  described  were  repeated  for  a  levee  cross  section 
with  irregular  geometry  using  a  standard  data  file  names  DATAIRR  (listed  in 
Appendix  D) .  The  screen  output  from  DATAIRR  using  default  parameters  is  shown 
in  Figure  7;  the  calculated  maximum  gradient  is  0.41.  In  this  figure,  the 
levee  geometry  display  has  been  windowed"  to  eliminate  some  of  the  pervious 
substratum  and  focus  on  the  area  near  the  levee.  As  evident  from  the  figure, 
this  file  models  a  top  blanket  having  a  thick  clay  plug  parallelling  the  land- 
side  toe,  and  a  broad,  water- filled  swale  on  the  landside.  Because  of  these 
irregularities,  no  closed  form  solution  can  be  obtained  for  this  problem. 
Results  of  the  parametric  study  are  shown  in  Figures  8  and  9.  It  is  seen  from 
Figure  8  that  node  spacings  of  10,  25,  and  50  ft  all  eventually  approach  a 
gradient  of  0.41  as  the  specified  tolerance  is  reduced;  however,  the  10  ft 
spacing  requires  reducing  the  tolerance  to  0.0001  ft  to  achieve  convergence. 

A  solution  of  this  problem  had  previously  been  obtained  (Wolff  1987)  using  the 
predecessor  program  LEVEEIRR.  The  earlier  solution  yielded  a  gradient  of 
0.37;  however,  LEVEEIRR  was  much  more  sensitive  to  the  way  nodes  were  speci¬ 
fied  and  did  not  ensure  close  node  spacing  near  the  levee  toe.  The  LEVEEMSU 
solution  is  considered  to  be  an  accurate  solution.  It  is  seen  from  Figure  9 
that  the  problem  DATAIRR  required  considerably  more  iterations  to  reach  a 
solution  than  did  DATACHK,  apparently  due  to  the  irregular  geometry.  The 
program  default  of  1,000  iterations  would  still  be  sufficient  to  achieve  con¬ 
vergence  with  the  default  tolerance  of  0.0005  ft  and  default  spacing  of  25  ft; 
however,  the  10  ft  node  spacing  discussed  above  required  over  4,000 
iterations . 
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Figure  9.  Number  of  iterations  versus  tolerance,  file  DATAIRR 


26.  Based  on  the  above  analyses,  program  default  parameters  were  set  at 
a  node  spacing  of  25  ft,  a  tolerance  of  0.0005  ft,  and  a  maximum  number  of 
iterations  of  1,000.  These  may  be  changed  by  the  user  during  program  execu¬ 
tion.  Finer  node  spacings  or  finer  tolerances  will  increase  the  number  of 
iterations;  the  user  is  cautioned  to  ensure  the  program  performs  sufficient 
iterations  to  reach  the  desired  solution. 

Comparison  to  Manual  Solution 

27.  The  program's  accuracy  was  checked  by  performing  a  manual  analysis 
for  the  conditions  modeled  by  the  input  file  DATACHK  (Appendix  C)  previously 
described.  The  program  output  is  shown  in  Appendix  C  and  the  manual  analysis 
is  shown  in  Appendix  E.  The  program  assumes  an  open  seepage  exit  at  the  last 
specified  section,  3,000  ft  landward  of  the  levee  toe.  In  the  manual  anal¬ 
ysis,  L3  distances  of  both  3,000  ft  and  infinity  were  checked.  The  computer 
solution  and  manual  solutions  are  compared  in  Figure  10,  which  plots  the  re¬ 
sidual  head  as  a  function  of  the  distance  from  the  lands ide  levee  toe  for  all 
three  solutions.  It  is  seen  that  the  computer  solution  is  accurate  and  that 
an  L3  distance  of  3,000  ft  accurately  models  an  infinitely  long  exit  condi¬ 
tion  for  this  case. 


Modeling  of  Finite  and  Infinite  Geometry 

28.  LEVEEMSU  always  models  open  entrance  and  exit  conditions  at  the 
first  and  last  specified  vertical  sections.  Infinitely  long  entrance  (L:) 
or  exit  (L3)  distances  must  be  approximated  by  specifying  the  beginning  or 
ending  sections  at  very  large  distances  from  the  levee.  To  investigate  how 
great  such  distances  should  be,  a  set  of  parametric  studies  was  performed  by 
systematically  altering  the  file  DATACHK  to  model  different  exit  lengths,  L3  . 
Results  are  shown  in  Figure  11,  in  which  the  gradient  at  the  landside  toe  is 
plotted  versus  L3  .  It  is  seen  that  results  become  constant  when  L3  ex¬ 
ceeds  2,000  to  3,000  ft,  or  in  this  case,  about  20  to  30  times  the  thickness 
of  the  pervious  substratum.  It  would  appear  that  such  a  ratio  should  accu¬ 
rately  model  infinitely  long  foundations;  however,  users  are  cautioned  to  make 
their  own  parametric  studies  for  cases  where  accuracy  is  critical.  Users  are 
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Figure  11.  Gradient  versus  Lo  ,  file  DATACHK 


is  a  math- 


further  cautioned  that  the  distance  to  effective  seepage  exit,  x3  , 
ematical  concept  used  in  conventional  analysis,  and  not  any  measure  of  levee 
geometry.  Specifying  the  last  vertical  section  at  some  calculated  x3  dis¬ 
tance  will  generally  result  in  calculated  gradients  that  are  too  low,  as  some 
seepage  always  in  fact  exits  the  blanket  beyond  the  x3  distance. 

Comparison  to  Program  LEVSEEP 

29.  A  computer  program  package  entitled  LEVSEEP  has  been  developed  under 
contract  to  W^S  by  Jaycor,  Inc.,  (Cunny,  Agostinelli,  and  Taylor  1989).  The 
program  LEVSEEP  performs  underseepage  analysis  for  the  traditional  model  of 
uniformly  thick  layers  and  uniform  properties.  The  program  also  performs  berm 
design  and  well  design  calculations  and  cost  estimates.  The  validation  for 
LEVSEEP  (Cunny,  Agostinelli,  and  Taylor  1989)  provided  comparative  solutions 
for  five  example  problems  worked  by  hand  and  by  computer.  Four  of  these  prob¬ 
lems  were  analyzed  using  LEVEEMSU  to  provide  further  program  verification. 

30.  Cunny 's  cross  section  No.  1  was  for  a  levee  on  a  pervious  foundation 
with  no  top  blanket.  As  LEVEEMSU  requires  that  a  top  blanket  be  present,  this 
cross  section  was  not  analyzed. 

31.  Cross  section  No.  2  was  for  a  levee  over  an  impervious  top  blanket 
and  a  foundation  with  open  entrance  and  exits  at  finite  distances  from  the 
levee.  This  cross  section  was  analyzed  using  LEVEEMSU  and  an  input  file  named 
JAYC0R2.  Results  are  shown  in  Figure  12.  As  LEVEEMSU  cannot  model  a  zero 
blanket  permeability,  the  blanket  was  modeled  with  a  permeability  of  1  x  10~7 
ft/min,  or  one  two-millionth  of  the  foundation  permeability.  As  modeled, 
LEVEEMSU  predicted  a  residual  head  of  8.392  ft  and  a  gradient  of  1.399;  the 
hand  analysis  predicted  a  head  of  8.333  ft  and  a  gradient  of  1.389.  The 
LEVEEMSU  results  exceeded  the  hand  analysis  by  0.7  percent.  Table  1  presents 
the  residual  head  and  gradient  calculated  by  LEVEEMSU,  LEVSEEP,  and  hand 
calculations . 

32.  Cross  section  No.  3  was  for  a  levee  over  a  semipervious  top  blanket 
having  a  foundation  with  a  finite  entrance  distance  and  an  infinite  exit  dis¬ 
tance.  This  cross  section  was  analyzed  using  an  input  file  named  JAYC0R3. 

The  finite  exit  distance  was  approximated  by  using  an  L3  distance  of  4,800 
ft.  R  u  ts  of  the  analysis  are  shown  in  Figure  13.  As  modeled,  LEVEEMSU 
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Table  1 

Comparison  of  LEVEEMSU.  LEVSEEP.  and  Hand  Calculations 
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predicted  a  residual  head  of  7.654  ft  and  a  gradient  of  1.276;  the  hand  analy¬ 
sis  predicted  a  head  of  7,723  ft  and  a  gradient  of  1.287.  The  LEVEEMSU  re¬ 
sults  were  lower  than  the  hand  analysis  by  0.9  percent. 

33.  Cross  section  No.  4  was  for  a  levee  over  a  semipervious  top  blanket 
having  a  foundation  with  a  finite  entrance  distance  and  exit  distances.  This 
cross  section  was  modeled  using  an  input  file  named  JAYC0R4.  Results  are 
shown  in  Figure  14.  As  modeled,  LEVEEMSU  predicted  a  residual  head  of  6.576 
ft  and  a  gradient  of  1.096.  The  hand  analysis  predicted  a  head  of  6.618  ft 
and  a  gradient  of  1.103.  The  LEVEEMSU  results  were  lower  than  the  hand  analy¬ 
sis  by  0.6  percent. 

34.  Cross  section  No.  5  was  for  Stovall,  MS,  Section  B-B.  This  section 
has  very  irregular  geometry  and  was  one  of  the  selected  prototype  test  sec¬ 
tions  for  LEVEEMSU  and  its  predecessor  LEVEEIRR.  The  hand  calculations  were 
made  for  Section  A-A  at  Stovall;  hence,  they  are  not  compared  in  Table  1. 
Analysis  of  this  section  is  shown  in  Part  IV  of  this  report. 

Effects  of  Blanket  Permeability 

35.  To  assess  the  consistency  of  program  behavior  with  respect  to  per¬ 
meability,  a  series  of  parametric  studies  were  performed.  The  data  file 
DATACHK  models  a  permeability  ratio  kf/kb  -  1,000  on  both  the  riverside  and 
landside  of  the  levee.  The  problem  was  altered  by  keeping  the  landside  per¬ 
meability  and  permeability  ratio  constant  (kf/kbl  -  1,000)  and  varying  the 
riverside  permeability  ratio  kf/kbr  from  1  to  1,000.  Then  the  riverside 
permeability  was  held  constant,  and  the  landside  permeability  ratio  changed  in 
a  similar  fashion.  Calculated  maximum  gradients  are  plotted  versus  permea¬ 
bility  ratio  in  Figure  15.  It  is  seen  that  the  program  exhibited  consistent 
and  expected  behavior,  with  the  gradient  increasing  with  increasing  riverside 
permeability  or  decreasing  landside  permeability,  and  vice  versa. 

Effects  of  Sloping  Ground 

36.  A  simple  but  perhaps  striking  illustration  of  the  effects  of  irregu¬ 
lar  geometry  is  the  case  of  a  sloping  ground  surface  landside  of  a  levee. 

This  is  a  commonly  encountered  condition  where  levees  are  founded  on  natural 
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Maximum  gradient  versus  permeability  ratio,  file  DATACHK 


levee  deposits.  To  assess  the  effects  of  ground  slope,  a  parametric  study  was 
made  by  successively  altering  the  file  DATACHK.  The  foundation  and  blanket 
thicknesses  were  kept  constant,  but  the  landside  ground  surface  slope  was 
varied  from  IV  on  1.000H  to  IV  on  50H  both  toward  the  levee  and  away  from  the 
levee.  An  example  screen  copy  from  the  study  is  shown  in  Figure  16.  Results 
are  shown  in  Figure  17.  The  observed  results  have  significant  implications. 

It  is  seen  that  a  ground  slope  as  little  as  h  to  1  percent  away  from  the  levee 
can  lower  the  gradient  several  tenths,  and  similar  slopes  toward  the  levee  can 
raise  the  gradient  a  similar  amount.  The  observed  differences  from  the  level 
ground  case  would  be  sufficient  to  result  in  sand  boils  at  sections  predicted 
to  be  safe,  and  safe  conditions  at  sections  predicted  to  have  boils.  Conven¬ 
tional  analysis  procedures  provide  no  means  to  assess  such  effects.  It  ap¬ 
pears  that  many  observed  inconsistencies  between  analysis  and  field  perform¬ 
ance  might  be  attributed  to  ground  slope  alone. 

Effects  of  Ditches 


37.  A  common  problem  in  underseepage  analysis  is  assessing  the  effects 
of  landside  ditches  and  determining  minimum  distances  that  ditches  should  be 
set  back  from  a  levee  to  provide  acceptable  gradients.  These  effects  can 
easily  be  assessed  using  LEVEEMSU.  To  illustrate  such  analyses,  a  file  named 
DATADCH  (Appendix  D)  was  created  and  systematically  modified  to  vary  the  dis¬ 
tance  between  the  landside  levee  toe  and  ditch  crown  and  to  vary  ditch  depth. 
In  all  cases,  the  ditch  had  a  10  ft  bottom  width  and  IV  on  3H  side  slopes. 

The  foundation  had  a  substratum  thickness  of  65  ft  and  a  top  stratum  thickness 
of  15  ft.  Both  constant  blanket  permeability  (PERMFLAG  =  "CONST")  and  vari¬ 
able  (PERMFLAG  =  "CURVE")  blanket  permeability  conditions  were  modeled.  A 
typical  screen  output  for  this  problem  is  shown  in  Figure  18.  Results  of  the 
study  are  summarized  in  Figure  19.  For  a  5-ft-deep  ditch  more  than  about  300 
to  350  ft  from  the  levee  toe,  the  gradient  at  the  levee  toe  (0.61)  exceeds  the 
gradient  at  the  ditch  for  the  conditions  modeled.  A  10-ft-deep  ditch  results 
in  excessively  high  gradients  even  at  distances  as  far  as  600  ft  from  the 
levees.  The  sections  modeled  with  variable  ("CURVE”)  blanket  permeabilities 
have  lower  permeability  values  except  at  the  ditch.  Using  the  curve  option, 
the  specified  permeability  is  for  a  10-ft-thick  blanket  in  the  ditch;  the 
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Maximum  gradient  versus  landside  ground  slope,  file  DATACHK 


Figure  18.  Copy  of  screen  output,  file  DATADCH 


5  ft  deep.  'CONST' 


Figure  19.  Maximum  gradient  versus  ditch  size  and  location 


program  adjusts  the  k  for  a  thickness  greater  than  10  ft  to  the  permeability 
for  a  10- ft  blanket  at  the  ditch  and  a  15- ft  blanket  away  from  the  ditch. 


Effects  of  Riverside  Borrow  Pits 


38.  Another  common  analysis  problem  is  assessing  the  effects  of  borrow 
pit  distance  and  depth  on  underseepage  conditions  for  new  levees  or  the  ef¬ 
fects  of  borrow  pit  enlargement  for  raising  existing  levees.  This  problem  is 
essentially  the  same  as  the  ditch  problem,  but  usually  alters  geometry  on  the 
riverside  of  the  levee.  Again,  LEVEEMSU  provides  a  tool  to  rationally  assess 
such  changes.  The  results  of  an  example  are  shown  in  Figure  20.  The  geometry 
of  file  DATACHK  was  modified  to  model  a  riverside  borrow  pit  300  ft  wide  and  5 
ft  deep  with  IV  on  3H  side  slopes  at  different  distances  from  the  levee  toe. 
These  results  are  shown  in  Figure  21.  It  is  seen  that  moving  a  borrow  pit 
closer  to  the  levee  increases  the  gradient  as  expected,  but  the  effect  is  much 
less  severe  than  cutting  ditches  on  the  landside. 

Effects  of  Relief  Wells 

39.  As  stated  in  Part  II,  LEVEEMSU  provides  the  capability  to  approxi¬ 
mately  model  the  effects  of  a  line  of  relief  wells  by  specifying  the  piezom¬ 
etric  head  at  one  landside  location.  The  program  calculates  the  well  flow  per 
foot  of  levee  required  to  reduce  the  average  piezometric  elevation  in  the  well 
line  to  the  specified  value.  The  designer  must  then  design  a  well  system 
consistent  with  these  results.  Relationships  obtained  between  piezometric 
elevations  and  well  flow  do  not  include  hydraulic  losses  or  partial  penetra¬ 
tion  effects  at  individual  wells.  They  should,  however,  be  useful  for  prelim¬ 
inary  assessments  of  the  need  for  wells  and  likely  numbers  and  spacing. 

40.  To  assess  the  program's  behavior,  two  parametric  studies  were  per¬ 
formed  using  the  input  file  DATAWELL,  which  is  essentially  the  file  DATACHK 
modified  to  specify  a  well  line  at  the  levee  toe.  In  the  first  study,  the 
specified  piezometric  elevation  in  the  well  line  was  varied.  Results  are 
shown  in  Figure  22.  In  this  figure,  both  maximum  gradient  and  well  flow  are 
plotted  versus  the  specified  average  head  at  the  well  line.  Results  of  such 
analysis  can  be  used  for  preliminary  design.  For  example,  if  it  is  desired  to 
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Figure  21.  Results  of  borrow  pit  study 


Gradient 


y/iud6  ‘mo|j  ||3M 


j  '^.uajpDJo  uunoiixD^ 


46 


Figure  22.  Maximum  gradient  and  well  flow  versus  average  head  in  well  line,  file  DATAWELL 


reduce  the  gradient  to  0.7,  the  figure  shows  that  a  well  system  must  be  de¬ 
signed  that  will  reduce  the  average  piezometric  elevation  the  well  line  to 
el  97.  Working  from  that  elevation,  the  system  must  be  capable  of  passing 
approximately  0.5  gal/min  per  ft  of  levee.  The  specified  piezometric  eleva¬ 
tion  resulting  in  zero  well  flow  is  98.87  ft,  matching  the  residual  head  of 
8.87  ft  obtained  for  the  analysis  of  DATACHK  without  wells.  In  the  second 
parametric  study,  the  specified  piezometric  elevation  was  maintained  constant 
and  the  foundation  permeability  was  varied.  Results  are  shown  in  Figure  23. 
As  expected,  well  flow  varies  in  proportion  to  the  foundation  permeability 
with  only  a  minor  change  in  gradient. 
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Figure  23.  Maximum  gradient  and  well  flow  versus  foundation  permeability,  file  DATAWELL 


PART  IV:  ANALYSIS  OF  PROTOTYPE  LEVEE  REACHES 


41.  Several  actual  levee  reaches  had  been  analyzed  using  LEVEEIRR,  the 
predecessor  of  LEVEEMSU,  to  provide  comparisons  of  predicted  versus  actual 
performance  (Wolff  1987).  Certain  of  these  reaches  have  been  reanalyzed  using 
LEVEEMSU.  These  reaches  had  previously  been  selected  as  having  relatively 
complete  and  reliable  piezometric  data  as  well  as  irregular  foundation 
conditions . 

42.  As  any  observed  piezometric  condition  can,  in  theory,  be  matched  to 
the  results  of  analysis  if  one  assumes  the  right  set  of  parameters,  these 
analyses  provide  only  a  partial  check  of  program  accuracy.  In  the  analyses 
reported  herein,  the  foundation  permeability  was  fixed  and  the  permeability 
values  for  the  landside  and  riverside  blankets  were  systematically  varied 
until  a  reasonable  match  was  obtained  between  predicted  and  measured  piezom¬ 
etric  data.  Assuming  the  program  provides  accurate  solutions,  these  analyses 
can  be  used  to  estimate  field  permeability  ratios,  landside  and  riverside, 
which  in  turn  can  be  used  in  the  design  of  seepage  control  measures. 

Rock  Island  District.  Hunt.  Piezometric  Range  B 

43.  This  piezometer  range  is  located  on  the  east  bank  of  the  Mississippi 
River  about  25  miles  upstream  of  Quincy,  IL,  in  the  pool  of  Lock  and  Dam 

No.  20.  Data  at  this  piezometer  range  has  previously  been  analyzed  by  Cunny 
(1980)  and  Wolff  (1987).  A  cross  section  of  the  site  is  shown  in  Figure  24; 
the  modeled  cross  section  is  shown  in  Figure  25.  Irregularities  in  the  pro¬ 
file  include  an  irregular  landside  ground  elevation  approximately  5  feet  high¬ 
er  than  the  riverside  and  a  blanket  of  variable  thickness. 

44.  Cunny 's  previous  analysis  assumed  uniformly  thick  blankets  (but 
different  thicknesses  and  different  permeabilities  on  opposite  sides  of  the 
levee) .  Cunny  found  that  observed  piezometric  conditions  at  the  levee  toe 
could  be  matched  using  permeability  ratios  of  k£/kbr  -  209  and  kf/kbl  -  64 
Thus,  kbl/kbr  -  3.26  .  The  previous  analysis  by  Wolff  (1987)  using  LEVEEIRR 
allowed  for  irregular  geometry  but  assumed  equal  blanket  permeabilities  on 
both  sides  of  the  levee.  This  analysis  showed  that  a  permeability  ratio 
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Riverside  Ltndside 


Modeled  cross  section,  Hunt 


kf/kb  between  20  and  64  provided  a  reasonable  match  to  observed  data,  both  at 
the  levee  toe  and  away  from  the  toe . 

45.  A  copy  of  a  typical  screen  output  from  the  LEVEEMSU  analysis  is 
shown  in  Figure  26.  The  foundation  permeability  was  taken  as  0.128  ft/min 
(640  x  10'4  cm  /sec).  A  riverside  permeability  ratio  of  kf/kbr  -  64  was  as¬ 
sumed  and  the  landside  permeability  ratio  was  varied.  Actual  piezometer  data 
represent  river  stages  of  492.54  ft  (1961),  496.13  ft  (1961),  and  499.0  ft 
(1973).  It  was  found  that  a  ratio  kf/kbl  —  4.27  provided  a  reasonable  match 
to  the  observed  conditions  at  piezometers  B-l  and  B-2,  as  shown  in  Figure  27. 
The  ratio  blanket  permeabilities,  kbl/kbr  is  found  to  be  15.  The  flat  re¬ 
sponse  of  piezometer  B-3  in  the  computer  model  arises  from  apparent  differ¬ 
ences  in  the  field  permeability  values  from  the  levee  toe  to  the  vicinity  of 
B-3.  The  blanket  must  be  modeled  as  very  pervious  to  match  observed  condi¬ 
tions  near  the  levee  toe;  when  this  is  done,  residual  heads  are  quite  low  at 
points  away  from  the  levee,  as  is  evident  from  Figure  26.  Reducing  the  land- 
side  blanket  permeability  would  improve  the  match  at  B-3  but  result  in  an 
overprediction  of  the  residual  head  near  the  levee  toe. 

46.  In  performing  this  analysis,  it  became  quite  apparent  that  predicted 
piezometric  elevations  remote  from  the  levee  are  quite  sensitive  to  the  speci¬ 
fied  landside  water  elevation.  Where  the  landside  water  elevation  is  modeled 
coincident  with  ridges  above  the  prevailing  ground,  the  model  will  produce 
downward  flow  from  the  ridges  to  the  aquifer,  causing  high  piezometric  levels. 
Where  the  ridged  are  modeled  as  being  above  the  landside  water  elevation,  more 
reasonable  results  are  obtained. 

47.  For  a  river  stage  at  the  crest  of  the  levee,  el  501.5  ft,  the  com¬ 
puter  predicts  a  maximum  gradient  of  0.385,  occurring  at  the  landside  levee 
toe . 

48.  Permeability  ratios  obtained  from  this  analysis  are  even  lower  than 
the  relatively  low  values  previously  calculated.  Since  a  perfect  match  be¬ 
tween  predicted  and  observed  piezometric  elevations  cannot  be  made  simultane¬ 
ously  for  all  piezometers  and  all  river  stages,  these  values  should  only  be 
considered  representative  of  the  actual  order  of  magnitude.  However,  one  can 
conclude  from  both  these  and  previous  analysis  that  the  field  ratio  kf/kbl 

is  in  the  range  of  4  to  80  and  the  ratio  kf/kbr  is  in  the  range  of  50  to 
250.  The  top  blanket  at  Hunt  appears  to  be  quite  pervious. 


52 


B-l  actual 


in 

c=> 

in 


in 


I 


f  +  t  +  * 

-  1 _ L_ 


in 

cn 


in 

OQ 


in 


(D 

> 

CC 


in 

co 


4 UOJ1-DA0 0IJ^8UU0Z8!d 


54 


Figure  27.  Piezometric  elevations  versus  river  elevation,  Hunt 


Memphis  District.  Commerce,  MS.  Line  H 

49.  This  piezometer  range  is  located  about  10  miles  north  of  Tunica,  MS. 
The  levee  is  located  about  2,200  ft  from  the  Mississippi  River  on  terrain 
characterized  by  numerous  ridges,  swales,  and  ditches.  Heavy  seepage  damage 
was  reported  during  the  1937  high  water.  Data  at  this  piezometer  range  has 
previously  been  analyzed  by  WES  (1964)  and  Wolff  (1987).  The  actual  and 
modeled  cross  sections  at  the  site  are  shown  superimposed  in  Figure  28.  Foun¬ 
dation  conditions  include  a  deep  riverside  borrow  pit  that  leaves  a  blanket  of 
only  a  few  feet  of  silt,  and  a  very  thick  clay  blanket  starting  about  1,600  ft 
landward  of  the  levee . 

50.  Previous  analyses  found  permeability  ratios  of  kf/kbr  <=  kf/kbl  =  580 
(WES  1964)  and  514  (Wolff  1987). 

51.  A  copy  of  a  typical  screen  output  from  LEVEEMSU  is  shown  in  Fig¬ 
ure  29.  The  foundation  permeability  was  taken  as  0.18  ft/min  (900  x  lO’4 
cm/sec) .  The  blanket  permeability  values  were  generated  using  the  variable 
permeability  option  described  in  Part  II  of  this  report.  On  the  riverside,  a 
curve  number  of  0.0007  was  specified,  corresponding  to  silt;  on  the  landside, 
a  curve  of  0.0005  was  specified,  corresponding  to  clay.  The  program  calcu¬ 
lates  the  blanket  permeability  at  each  node  using  the  blanket  thickness  at 
that  node.  Results  of  the  present  analysis  are  compared  with  actual  piezom¬ 
etric  data  in  Figure  30;  piezometer  data  correspond  to  river  stages  of  202.7 
and  205.0  ft  (May  18  and  22,  1961).  It  is  of  interest  to  note  that  the  vari¬ 
able  permeability  option  provided  a  reasonably  good  match  to  actual  perform¬ 
ance  on  the  first  try. 

52.  For  a  river  stage  at  the  levee  crest,  el  220.2  ft,  the  computer 
model  predicts  a  maximum  gradient  of  1.06  at  x  =  330  ft  landside  of  the 
centerline,  the  location  of  the  landside  berm  toe  and  piezometer  9-x. 

Memphis  District.  Stovall.  MS.  Piezometric  Line  B 

53.  This  piezometer  range  is  located  about  3.5  miles  west  of  Stovall, 

MS.  Seepage  damage  occurred  at  the  site  during  the  1937  high  water.  Founda¬ 
tion  conditions  are  shown  superimposed  on  the  modeled  conditions  in  Figure  31. 
These  conditions  are  an  example  of  extremely  dissimilar  landside  and  riverside 
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DISTANCE  FROM  LEVEE  IN  FEET 


Actual  and  modeled  cross  section,  Commerce  (after  WES,  1964) 
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Figure  30.  Piezometric  elevations  versus  river  elevation,  Commerce 


DISTANCE  FROM  LEVEE  IN  FEET 


Actual  and  modeled  cross  section,  Stovall  (after  WES,  1964) 


blanket  conditions.  The  riverside  blanket  is  predominantly  silt  and  much  of 
it  has  been  excavated  for  borrow.  The  riverside  blanket  is  clay,  typically  15 
ft  thick.  A  seepage  berm  has  been  constructed  on  the  landside.  Piezometric 
data  are  available  for  May  1961  (WES  1964) . 

54.  Previous  analysis  by  the  WES  (1964)  obtained  a  ratio  kf/kbl  of  600. 
Previous  analysis  using  LEVEEIRR  (Wolff  1987)  indicated  kf/kbl  ratios  in  the 
range  of  432  to  1,000;  however,  the  modeling  of  equal  blanket  permeabilities 
on  both  sides  of  the  levee  did  not  provide  a  match  to  observed  data  at  several 
piezometers  simultaneously. 

55.  Using  LEVEEMSU,  the  section  was  analyzed  using  the  variable  permea¬ 
bility  option  described  in  Part  II.  The  best  match  to  observed  conditions  was 
obtained  using  a  riverside  curve  number  of  0.0006  corresponding  to  a  silty 
clay  or  clayey  silt,  and  a  landside  curve  number  of  0.0004,  corresponding  to  a 
clay . 

56.  A  copy  of  a  screen  display  of  the  modeled  cross  section  is  shown  in 
Figure  32.  Note  that  landside  and  riverside  are  reversed  from  the  previous 
figure.  The  entire  riverside  geometry  extends  more  than  a  mile  to  the  river; 
however,  the  display  has  been  "windowed"  to  show  more  detail  near  the  levee. 
The  permeability  curve  numbers  stated  above  were  found  to  yield  predicted 
piezometric  conditions  that  match  observed  data  from  four  piezometers  at 
greatly  different  distances  from  the  levee,  generally  within  less  than  1  ft. 
Furthermore,  the  match  was  obtained  with  only  a  few  trials,  by  making  minor 
adjustments  in  the  curve  number.  The  results  of  the  computer  model  are  com¬ 
pared  to  observed  conditions  in  Figure  33.  From  this  analysis,  it  appears 
that  the  proposed  approximations  of  the  LMVD  permeability  versus  blanket 
thickness  curves  yield  reasonable  solutions  when  used  with  reasonably  accurate 
geometry  data. 

Vicksburg  District.  Bolivar.  MS.  Piezometric  Line  D 

57.  This  piezometer  range  is  located  along  the  east  bank  levee  of  the 
Mississippi  River  2  miles  north  of  Benoit,  MS.  The  river  at  this  site  is 
about  8  miles  from  the  levee;  however,  Bolivar  Chute  lies  about  1,200  to  1,500 
ft  riverward  of  the  levee.  A  range  of  nine  piezometers,  D-l  through  D-9,  runs 
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INPUT  FILE:  stovall -  210  - 05-31-1989 

Memphis  District  —  Stovall  13:28:35 
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Figure  32.  Copy  of  screen  output,  Stovall 


11  actual 


Piezometric  elevations  versus  river  elevation,  Stovall 


perpendicular  to  the  levee  at  this  site.  A  cross  section  of  foundation  condi¬ 
tions  with  the  modeled  cross  section  superimposed  is  shown  in  Figure  34.  Ir¬ 
regularities  in  the  profile  include  riverside  borrow  pits,  landside  sublevees, 
a  landside  ditch,  and  a  massive  clay- filled  abandoned  channel  beginning  about 
1,000  ft  landside  of  the  levee. 

58.  Previous  analysis  by  the  WES  (1964)  found  kf/kbl  ratios  in  the 
range  of  100  to  200.  Previous  analyses  by  Wolff  (1987)  found  that  a  ratio 
kf/kbi  of  1,000  best  fit  the  data;  however,  little  difference  was  noted  for  a 
value  of  100. 

59.  For  the  present  analysis,  a  foundation  permeability  of  0.24  ft/min 
(1,200  x  10"*  cm/sec)  was  assumed.  A  copy  of  a  typical  screen  output  from 
LEVEEMSU  is  shown  in  Figure  35.  Actual  piezometer  data  represent  river  stages 
of  147.0  ft  (1961)  and  151.7  ft  (1973).  More  weight  was  given  to  the  1973 
data  in  performing  the  analysis.  It  was  found  that  ratios  of  k£/kbl  =  267 
and  kf/kbr  —  1,000  provided  a  reasonable  match  to  observed  conditions  as 
shown  in  Figures  36  and  37.  The  ratio  of  blanket  permeabilities,  kbl/kbr  is 
found  to  be  3.75. 

60.  For  a  river  stage  at  the  project  flow  line  of  el  166.4  ft,  the  com¬ 
puter  model  predicts  a  maximum  gradient  of  0.66,  occurring  at  the  levee  toe. 

61.  Results  of  the  study  again  support  the  observation  that  landside 
blanket  permeability  values  during  high  water  are  several  times  higher  than 
riverside  values. 
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Figure  34.  Actual  and  modeled  cross  section,  Bolivar 
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Figure  36.  Piezometric  elevations  versus  river  elevation,  Bolivar 


PART  V:  DESIGN  APPLICATIONS 


62.  Although  the  program  LEVEEMSU  was  developed  primarily  as  an  analysis 
tool,  it  can  be  employed  quite  readily  for  preliminary  design  of  underseepage 
control  measures  such  as  seepage  berms  and  relief  wells. 

Design  of  Seepage  Berms 

63.  If  a  semipervious  berm  is  assumed  to  have  a  vertical  permeability  or 
permeability  curve  number  equal  to  that  of  the  landside  top  blanket,  a  berm 
can  be  modeled  by  adjusting  the  geometry  of  the  landside  blanket  to  include 
the  berm.  By  a  few  trial  and  error  adjustments,  a  designer  can  size  a  berm 
that  will  result  in  any  desired  gradient  at  any  landside  location.  Figure  38 
illustrates  a  copy  of  the  screen  display  for  an  input  file  entitled  DATABERM. 
This  file  was  developed  by  adjusting  the  geometry  of  file  DATACHK  to  include  a 
300-ft-long  berm  with  a  thickness  of  5  ft  at  the  levee  toe  and  2  ft  at  the 
berm  toe.  The  maximum  gradient  of  0.71  occurs  at  the  berm  toe  (x  =  2,106). 
From  the  printed  output  file,  the  gradient  at  the  levee  toe  and  any  other  lo¬ 
cation  can  be  obtained.  Furthermore,  the  graphic  display  of  the  relationship 
of  the  piezometric  line  to  the  berm  surface  provides  the  designer  a  clear 
basis  for  making  trial  and  error  adjustments  of  berm  geometry. 

64.  A  logical  improvement  to  the  program  would  be  to  automate  the  berm 
design  process.  The  program  could  be  given  a  set  of  allowable  gradient  values 
for  different  distances  from  the  levee  and  would  automatically  adjust  the  berm 
surface  until  the  desired  gradients  were  attained. 

Design  of  Relief  Wells 

65.  As  previously  demonstrated  in  Part  III,  use  of  the  option  to  specify 
piezometric  head  at  one  node  allows  approximate  modeling  of  a  well  line  and 
provides  the  designer  with  a  relationship  between  average  piezometric  eleva¬ 
tion  at  the  well  line  and  well  flow  per  foot  of  levee.  These  results  may  be 
sufficient  for  preliminary  studies  and  cost  estimates  where  only  the  approxi¬ 
mate  number  of  wells  is  required.  Detailed  design  of  a  well  system  including 
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spacing,  hydraulic  losses,  and  partial  penetration  effects  is  beyond  the  scope 
of  the  present  program. 
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PART  VI:  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


Summary 

66.  A  software  package  for  underseepage  analysis  of  levees,  LEVEEMSU, 
has  been  developed  and  tested.  The  program  allows  an  explicit  description  of 
the  actual  geometry  of  the  top  blanket,  pervious  substratum,  and  landside 
water,  which  need  not  be  regular  or  uniform.  Blanket  permeability  may  be  con¬ 
stant  or  a  function  of  blanket  thickness.  The  residual  head  and  gradient  are 
calculated  along  the  entire  length  of  the  top  blanket  using  a  finite  different 
approximation  of  Bennett's  equation.  A  graphic  display  of  the  problem  and 
results  are  provided. 

67.  A  number  of  parametric  studies  were  performed  to  test  the  program 
and  demonstrate  it  capabilities,  which  include  analysis  of  levee  sections  hav¬ 
ing  ditches,  borrow  pits,  sloping  ground,  and  other  irregularities.  Several 
prototype  reaches  were  analyzed  to  assess  and  demonstrate  the  program's  capa¬ 
bility  to  match  observed  piezometric  conditions  along  an  entire  cross  section 
or  piezometer  range. 


Conclusions 


68.  Although  the  capabilities  of  LEVEEMSU  would  be  inherent  in  a 
general-purpose  finite  element  program,  the  development  and  use  of  this 
special-purpose  program  should  have  special  advantages  for  levee  analysis  and 
design: 

a.  The  program  is  relatively  short  and  can  be  run  on  virtually  any  IBM 
compatible  microcomputer  running  the  MS  DOS  (TM)  operating  system 
and  having  CGA  or  EGA  graphics  capabilities. 

b.  The  program  input  and  output  is  specifically  oriented  to  levees. 

For  example,  information  on  the  gradient  through  the  top  blanket 
would  have  to  be  manually  extracted  from  the  output  of  a  general- 
purpose  program. 

o.  The  assumptions  made  in  the  analysis  are  identical  to  the  assump¬ 
tions  inherent  in  conventional  analysis.  Thus,  program  solutions 
should  allow  a  user  to  match  conventional  analyses  and  then  extend 
them  to  more  complex  actual  conditions. 

69.  Based  on  extensive  program  testing  and  parametric  studies,  a  number 
of  conclusions  are  drawn: 
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a.  The  program's  algorithm  for  automatic  generation  of  nodes  at  vari¬ 
able  spacing  provides  considerably  more  accurate  and  consistent  so¬ 
lutions  than  the  predecessor  program,  LEVEEIRR. 

b.  The  program  solutions  match  hand  solutions  for  cases  of  uniform  ge¬ 
ometry  where  companion  solutions  can  be  obtained. 

c.  Parametric  studies  show  that  the  program  exhibits  consistent  and 
reasonable  behavior  with  respect  to  changes  in  permeability,  ground 
slope,  ditch  location,  relief  well  characteristics,  and  other 
variables . 

d.  Results  of  the  ground  slope  study  suggest  that  many  observed  dis¬ 
crepancies  between  actual  and  measured  gradients  might  be  attribut¬ 
able  to  ground  slope  effects  alone. 

e.  Results  of  the  prototype  reach  analyses  show  that  field  permeability 
ratios  can  be  estimated  by  systematically  varying  program  input  un¬ 
til  a  reasonable  match  to  observed  conditions  is  obtained.  Using 
the  program,  this  can  be  done  more  precisely  than  by  conventional 
analysis  as  the  analyst  does  not  need  to  assign  constant  "design" 
values  to  parameters  and  dimensions  that  are  in  fact  variable. 

f .  LEVEEMSU  provides  a  useful  and  convenient  analysis  and  design  tool 
that  should  allow  designers  to  more  accurately  model  actual  condi¬ 
tions.  Presumably,  this  should  lead  to  more  accurate  predictions. 

70.  Flood  protection  is  a  complex  system  involving  design,  construction, 
maintenance,  and  performance  evaluation  of  levees.  An  analysis  program  such 
as  LEVEEMSU  can  be  but  one  link  in  such  a  system.  If  experience  proves  that 
the  program  provides  a  capability  to  more  accurately  evaluate  and  predict  un¬ 
derseepage  conditions,  this  should  lead  to  a  reevaluation  of  design  crit-eria 
with  a  view  to  both  reducing  cost  and  improving  safety. 


Recommendations 


71.  Based  on  the  results  of  this  research,  the  following  recommendations 
are  made : 

a.  The  program  LEVEEMSU  should  be  field  tested  by  use  in  District  of¬ 
fices  and  the  need  for  any  corrections  or  improvements  assessed. 

b.  The  improvement  of  LEVEEMSU  to  provide  a  design  mode  should  be  con¬ 
sidered.  Given  a  set  of  berm  design  criteria,  the  program  could 
adjust  berm  dimensions  until  the  desired  gradients  are  obtained. 

c.  Development  of  an  extended  version  of  LEVEEMSU  to  analyze  three - 
layer  foundations  should  be  investigated.  In  many  instances,  a  lay¬ 
er  of  silty  sand  or  sandy  silt  between  the  substratum  and  top  blan¬ 
ket  may  have  significant  vertical  and  horizontal  flow  components. 

In  such  sections,  construction  of  a  reasonable  analysis  model  based 
on  two  layers  is  known  to  be  difficult. 
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d.  The  need  for  refining  levee  design  criteria  should  be  assessed. 

Many  current  criteria,  such  as  dimensions  and  location  of  borrow 
pits  and  ditches  are  necessarily  arbitrary  and  conservative  due  to 
the  lack  of  a  rational  analysis  procedure.  LEVEEMSU  provides  the 
capability  for  site-specific  analysis  of  such  items.  Likewise,  cri¬ 
teria  for  sizing  seepage  berms  presumably  includes  an  allowance  for 
discrepancies  between  the  actual  ground  conditions  and  the  ground 
conditions  analyzed.  The  use  of  LEVEEMSU  reduces  such 
discrepancies . 
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APPENDIX  A:  DATA  INPUT  FOR  LEVEEMSU 


1.  Input  data  are  read  from  a  standard  ASCII  text  file  which  can  be  ere 
ated  using  a  word  processor  or  text  editor.  Information  in  the  input  file 
includes  coordinates  of  the  soil  profile  and  information  on  material  proper¬ 
ties.  Coordinate  information  is  specified  in  terms  of  vertical  sections  cut 
through  the  profile.  At  each  specified  vertical  section,  the  base  of  the  per 
vious  substratum,  base  of  the  top  blanket,  and  ground  surface  are  specified. 
For  landside  sections,  the  landside  water  surface  is  also  specified.  Layer 
boundaries  are  assumed  to  vary  linearly  between  specified  sections. 

2.  An  example  input  file  with  corresponding  variable  names  and  defini¬ 
tions  is  shown  below;  a  sketch  of  the  corresponding  cross  section  is  shown  in 
Figure  A-l. 


Example  Data  File: 

IRREGULAR  FOUNDATION 

TEST  PROBLEM 

0.200 

2  "CONST"  .0002  175 
750.0  60.0  140.0  158.0 
1750.0  60.0  140.0  160.0 
4  "CONST"  .0002 
1900.0  60.0  140.0  160.0  160.0 
2400.0  60.0  120.0  158.33  158.33 
2800.0  62.0  140.0  155.0  158.33 
4900.0  70.0  140.0  158.0  158.33 
NO  WELLS 


Variable  Names 

TITLE1$ 

TITLE2$ 

KF 

NRIVSECS  PERMFLAGR$  PERMRIV  YRIV 
X(l) ,  Yl(l) ,  Y2(l) ,  Y3 (1) 

II 

NLANDSECS  PERMFLAGL  PERMLAND 
X(*),  Y1 (*) ,  Y2  (*) ,  Y3(*) ,  YWATER(*> 

If 

It 

II 

WELLFLAG$ 


(If  it  is  desired  to  model  the  effect  of  a  line  of  relief  well  line,  the  last 
line  above  is  deleted  and  two  lines  are  added  as  shown  below) 


WELL  WELLFLAG$ 

1920  162  XWELL  YWELL 


A-l 


Definition  of  input  variables 


Definition  of  Variables 


3.  Variable  definitions  are  given  below  in  the  order  they  are  listed  in 
the  sample  input  file  above. 

TITLE1$  is  the  first  title  line. 

TITLE2$  is  the  second  title  line  (two  are  required). 

KF  is  the  horizontal  permeability  of  the  pervious  substratum,  in  ft/min. 

NRIVSECS  is  the  number  of  vertical  sections  used  to  describe  the  problem 
geometry  on  the  riverside  of  the  levee.  NRIVSECS  must  be  two  (2)  or 
more.  The  first  section  is  an  open  entrance;  the  last  section  is  the 
riverside  toe  of  the  levee. 

PERMFLAGR$  is  a  flag  that  indicates  how  the  riverside  permeability  is  to 
be  specified.  Use  the  value  "CONST"  or  "const"  to  specify  a  constant 
riverside  blanket  permeability.  Use  the  value  "CURVE"  or  "curve"  to  cal¬ 
culate  the  riverside  blanket  permeability  as  a  function  of  the  blanket 
thickness  z  . 

PERMRIV:  If  PERMFLAGR$  is  "CONST”,  PERMRIV  is  the  vertical  permeability 

of  the  riverside  top  blanket  in  ft/min.  If  PERMFLAG$  is  "CURVE",  PERMRIV 
is  the  vertical  permeability  for  a  blanket  thickness  of  10  ft,  and  the 
program  will  calculate  the  permeability  for  other  thicknesses  using  the 
method  described  in  Part  II. 

YRIV  is  riverside  water  elevation. 

X(l) ,  Yl(l),  Y2(l),  Y3(l)  are  the  geometry  data  for  the  first  vertical 
section,  with  x  increasing  from  riverside  to  landside. 

X(l)  is  the  x  coordinate 

Yl(l)  is  the  base  of  the  pervious  substratum 

Y2(l)  is  the  top  of  the  pervious  substratum/base  of  the  top  blanket 
Y3(l)  is  the  top  of  the  ground 

These  lines  are  repeated  for  each  riverside  vertical  section.  (NRIVSECS  lines 
in  all) . 

NLANDSECS  is  the  number  of  vertical  sections  used  to  describe  the  problem 
geometry  on  the  landside  of  the  levee.  NLANDSECS  must  be  two  (2)  or 
more.  The  first  section  is  the  landside  toe  of  the  levee;  the  last  sec¬ 
tion  is  an  open  exit. 

PERMFLAGR$  is  a  flag  that  indicates  how  the  landside  permeability  is  to 
be  specified.  Use  the  value  "CONST"  or  "const"  to  specify  a  constant 
landside  blanket  permeability.  Use  the  value  "CURVE"  or  "curve"  to  cal¬ 
culate  the  landside  blanket  permeability  as  a  function  of  the  blanket 
thickness  z  . 

PERMLAND:  If  PERMFLAGR$  is  "CONST",  PERMLAND  is  the  vertical  permeabil¬ 

ity  of  the  landside  top  blanket  in  ft/min.  If  PERMFLAG$  is  "CURVE”, 
PERMLAND  is  the  vertical  permeability  for  a  blanket  thickness  of  10  ft, 


and  the  program  will  calculate  the  permeability  for  other  thicknesses 
using  the  method  described  in  Part  II. 

X(*) ,  Yl(*) ,  Y2(*),  Y3(*),  YWATER(*)  are  the  geometry  data  for  the  first 
landside  vertical  section  (at  the  landside  toe). 

X(*)  is  the  x  coordinate 

Yl(*)  is  the  base  of  the  pervious  substratum 

Y2(*)  is  the  top  of  the  pervious  substratum/base  of  the  top  blanket 
Y3(*)  is  the  top  of  the  ground 

YWATER(*)  is  the  elevation  of  the  free  water  surface  (typically) 
equal  to  or  above  the  ground  surface) . 

These  lines  are  repeated  for  each  landside  vertical  section.  (NLANDSECS  lines 
in  all) . 

WELLFLAG$  is  a  flag  which  tells  the  program  to  read  one  additional  line 
giving  a  specified  piezometric  elevation  at  one  location,  which  can  be 
used  to  simulate  a  line  of  relief  wells.  If  a  relief  well  (or  specified 
piezometric  head)  is  to  be  specified,  enter  the  word  WELL  or  well  on  this 
line.  If  this  option  is  not  desired,  enter  any  other  word(s) ,  such  as  NO 
WELLS,  STOP  or  END. 

The  following  variables  are  only  used  if  WELLFLAG$  is  WELL  or  well: 

XWELL  is  the  x  coordinate  where  the  piezometric  elevation  is  to  be 
specified.  If  a  node  is  not  generated  at  this  location,  the  program  will 
m<~ve  it  to  the  nearest  node.  It  is  recommended  that  this  value  be  the 
same  as  the  x  coordinate  of  one  of  the  specified  landside  sections. 

YWELL  is  the  y  coordinate  of  the  specified  piezometric  elevation  at 
XWELL.  The  piezometric  elevation  will  be  forced  to  the  specified  value 
at  the  specified  location.  It  is  equivalent  to  the  average  piezometric 
elevation  in  a  line  of  wells . 
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APPENDIX  B:  RUNNING  LEVEEMSU 


1.  LEVEEMSU.EXE  is  a  stand-alone  executable  file  compiled  from  a  high- 
level  BASIC  source  code.  The  program  runs  on  IBM  compatible  personal  compu¬ 
ters  using  the  MS  DOS  operating  system  and  having  CGA  or  EGA  graphics  capabil¬ 
ities.  A  math  coprocessor  is  not  required,  but  the  program  runs  considerably 
faster  if  a  coprocessor  is  installed.  A  graphics  printer  is  required  to  ob¬ 
tain  a  printer  plot  of  the  foundation  geometry  and  results. 

2.  Before  running  the  program,  prepare  and  save  one  or  more  data  files 
in  standard  ASCII  format  using  a  word  processor  or  text  editor.  Format  for 
data  files  is  given  in  Appendix  A.  Any  number  of  files  can  be  analyzed  with¬ 
out  exiting  the  program. 

3.  To  obtain  a  printer  plot  of  a  CGA  screen  display,  the  DOS  command 
GRAPHICS  (program  GRAPHICS ,. COM)  must  be  resident  on  the  system  and  executed 
before  running  the  program.  For  example  (User  input  is  underlined) : 

C>  GRAPHICS 


To  obtain  a  printer  plot  of  an  EGA  screen  display  from  monochrome  mode, 
the  public  domain  program  EPSON.COM  or  a  similar  EGA  screen  dump  program 
should  be  executed  before  running  the  program.  For  example: 


C>  EPSON 

4.  To  run  the  program,  log  to  the  drive  where  the  program  resides  (usu¬ 
ally  drive  C:)  and  type  the  program  name: 

C>  LEVEEMSU 


The  program  will  display  introductory  information  and  ask  what  type  of 
graphics  display  is  available.  Respond  with: 


EGAC  or  EGAM  or  CGA 


The  option  EGAC  will  provide  a  high  resolution  color  graphic  display  of 
the  problem.  The  option  EGAM  will  provide  a  high  resolution  monochrome 
display,  suitable  for  copying  to  a  graphics  printer.  The  option  CGA  will 
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provide  a  medium  resolution  monochrome  graphic  display  of  the  problem, 
suitable  for  copying  to  a  graphics  printer. 


5.  The  program  will  ask  for  the  input  file  name.  Enter  the  file  name. 
Prefix  the  file  name  with  the  drive  identifier  and/or  path,  if  different  from 
the  default  drive  and  directory.  Include  the  extension  if  present.  For 
example : 

DATAIRR 

A: DATAIRR 

A: DATAIRR .DAT 

are  all  valid  data  file  identifiers.  If  the  input  file  specified  does  not 
exist,  the  program  will  prompt  for  the  file  name  again. 

6.  The  program  will  ask  for  the  output  file  name  in  a  similar  fashion. 

If  the  named  file  does  not  exist,  it  will  be  created.  If  it  does  exist,  it 
will  be  overwritten. 

7.  The  program  will  then  ask  if  you  wish  to  change  any  of  the  default 
settings  for  closure  tolerance,  maximum  iterations,  and  maximum  node  spacing 
near  the  levee  toe.  If  you  are  satisfied  with  the  present  settings,  press 
return  for  each  value.  If  you  wish  to  change  these  values,  type  in  the  new 
value , 

8.  The  program  will  print  a  summary  of  your  input  data.  If  you  wish  a 
hard  copy,  press  Shift-PrtSc  before  pressing  Return  to  continue.  Input 
data  are  also  saved  to  a  specified  output  file, 

9.  The  program  will  provide  a  graphic  display  of  your  input  data.  If 
you  wish  a  hard  copy,  press  Shift-PrtSc  before  pressing  Return  to 
continue . 

10.  The  program  will  then  solve  for  the  head  and  gradient  along  the  base 
of  the  blanket.  When  the  solution  is  complete,  it  will  plot  the  piezometric 
grade  line.  Ii  you  wish  a  hard  copy,  press  Shift-PrtSc  before  pressing 
Return  to  continue. 


R-2 


11.  The  program  will  ask  whether  the  graphics  window  is  to  be  changed. 
This  feature  allows  zooming  in  on  a  particular  area  of  interest.  Enter  Y(es) 
or  N(o)  .  If  the  response  is  yes,  the  program  will  display  the  coordinates 
of  the  current  window  and  prompt  you  for  the  new  window  coordinates.  These 
are  entered  as  minimum  and  maximum  x  and  y  values,  respectively,  sepa¬ 
rated  by  commas .  The  window  boundaries  can  be  changed  as  often  as  desired. 

12.  The  program  will  ask  whether  the  results  are  to  be  listed  to  the 
printer.  Enter  Y(es)  or  N(o)  .  Whether  printed  or  not,  results  are  saved 
in  the  output  file  for  later  printing  using  a  word  processor  or  the  DOS  COPY 
or  PRINT  commands . 

13.  The  program  will  then  prompt  for  a  new  problem.  Enter  Y(es)  to  go 
to  a  new  problem,  or  N(o)  to  quit. 
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APPENDIX  C :  EXAMPLE  RUN ,  FILE  DATACHK 


-  LEVEEMSU  - 

UNDERSEEPAGE  ANALYSIS 
FOR 

IRREGULAR  FOUNDATION  CONDITIONS 

Thomas  F.  Wolff 
Michigan  State  University 


U.  S.  Army  Engineer 

Waterways  Experiment  Station  | 

Geotechnical  Laboratory 

I 


Last  revision  May  31,  1989 


Enter  graphics  mode 
EGAC  for  Hi  Res  color 
EGAM  for  Hi  Res  monochrome 
CGA  for  Med  Res  monochrome 
?  egac 


ENTER  INPUT  FILE  NAME 
?  datachk 

ENTER  OUTPUT  FILE  NAME 
?  out 
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DEFAULT  SETTINGS 

Press  return  if  OK 

Type  new  number  to  change 

Tolerance  =  .0005 

7 

Maximum  iterations  =  1000 

r> 

Maximum  node  spacing  near  levee  = 


Data  file:  datachk 

IRREGULAR  FOUNDATION 
TEST  PROBLEM 

Foundation  Fermeability 

.2 


Riverside 

sections 

permf 1 ag 

perm 

2 

CONST 

.0002 

X 

0 

l1 

ll 

1500 

0 

80 

Landside 

sections 

2 

permf lag 
CONST 

perm 

.0002 

X 

V1 

y2 

1800 

0 

80 

5000 

0 

80 

Press  Return  to  Continue 


yri  v 
110 


90 


ywater 

90 

90 
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Listed  Output  File 


PROGRAM  LEVEEHSU  Nay  31,  198?  edition 
INPUT  FILE:  datachk 
OUTPUT  FILE:  out 

IRREGULAR  FOUNDATION 
TEST  PROBLEM 


KF  =  .2 

PERNFLAGR  =  CONST  PERHRIV  -  .0002 
PERNFLAGL  =  CONST  PERNLAND  =  .0002 


XX 

yyl 

yy2 

yy3 

yyuater 

d 

z  kb 

0.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

500.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

700.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

900.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1000.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1100.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1200.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1250.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1300.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1350.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1400.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1425.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1450.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1475.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1500.00 

0.00 

80.00 

90.00 

110.00 

80.0 

10.0  0.20000E-03 

1525.00 

0.00 

80.00 

90.00 

108.33 

80.0 

10.0  O.OOOOOE+OO 

1550.00 

0.00 

80.00 

90.00 

106.67 

80.0 

io.o  o.oooooe+oo 

1575.00 

0.00 

80.00 

90.00 

105.00 

80.0 

10.0  O.OOOOOE+OO 

1600.00 

0.00 

80.00 

90.00 

103.33 

80.0 

10.0  0. OOOOOE+OO 

1625.00 

0.00 

80.00 

90,00 

101.67 

80.0 

10.0  O.OOOOOE+OO 

1650.00 

0.00 

80.00 

90.00 

100.00 

80.0 

10.0  O.OOOOOE+OO 

1675.00 

0.00 

80.00 

90.00 

98.33 

80.0 

10.0  O.OOOOOE+OO 

1700.00 

0.00 

80.00 

90.00 

96.67 

80.0 

10.0  O.OOOOOE+OO 

1725.00 

0.00 

80.00 

90.00 

95.00 

80.0 

10.0  O.OOOOOE+OO 

1750.00 

0.00 

80.00 

90.00 

93.33 

80.0 

10.0  O.OOOOOE+OO 

1775.00 

0.00 

80.00 

90.00 

91.67 

80.0 

10.0  O.OOOOOE+OO 

1800.00 

0.00 

80.00 

90.00 

90.00 

80.0 

10.0  0. 20000E-03 

1825.00 

0.00 

80.00 

90.00 

90.00 

80.0 

10.0  0.20000E-03 

1850.00 

0.00 

80.00 

90.00 

90.00 

80.0 

10.0  0.20000E-03 

1875.00 

0.00 

80.00 

90.00 

90.00 

80.0 

10.0  0.20000E-03 

1900.00 

0.00 

80.00 

90.00 

90.00 

80.0 

10.0  0.20000E-03 

1925.00 

0.00 

80.00 

90.00 

90.00 

80.0 

10.0  0.20000E-03 

1975.00 

0.00 

80.00 

90.00 

90.00 

80.0 

10.0  0.20000E-03 

2025.00 

0.00 

80.00 

90.00 

90.00 

80.0 

10.0  0.20000E-03 

2075.00 

0.00 

80.00 

90.00 

90.00 

80.0 

10.0  0.20000E-03 
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2125.00 

0.00 

80.00 

90.00 

90.00 

80.0 

2225.00 

0.00 

80.00 

90.00 

90.00 

80.0 

2325.00 

0.00 

80.00 

90.00 

90.00 

80.0 

2425.00 

0.00 

80.00 

90.00 

90.00 

80.0 

2625.00 

0.00 

80.00 

90.00 

90.00 

80.0 

2825.00 

0.00 

80.00 

90.00 

90.00 

80.0 

3325.00 

0.00 

80.00 

90.00 

90.00 

80.0 

3025.00 

0.00 

80.00 

90.00 

90.00 

80.0 

4325.00 

0.00 

80.00 

90.00 

90.00 

80.0 

4825.00 

0.00 

80.00 

90.00 

90.00 

80.0 

5000.00 

0.00 

80.00 

90.00 

90.00 

80.0 

XX 

piesel 

reshead 

z 

i 

0.00 

110.00 

0.00 

10.00 

0.000 

500.00 

108.33 

-1.67 

10.00 

-0.167 

700.00 

107.41 

-2.59 

10.00 

-0.259 

900.00 

106.36 

-3.64 

10.00 

-0.364 

1000.00 

105.77 

-4.23 

10.00 

-0.423 

1100.00 

105.12 

-4.88 

10.00 

-0.488 

1200.00 

104.42 

-5.58 

10.00 

-0.558 

1250.00 

104.04 

-5.96 

10.00 

-0.596 

1300.00 

103.64 

-6.36 

10.00 

-0.636 

1350.00 

103.23 

-6.77 

10.00 

-0.677 

1400.00 

102.79 

-7.21 

10.00 

-0.721 

1425.00 

102.56 

-7.44 

10.00 

-0.744 

1450.00 

102.33 

-7.67 

10.00 

-0.767 

1475.00 

102.09 

-7.91 

10.00 

-0.791 

1500.00 

101.84 

-8.16 

10.00 

-0.816 

1525.00 

101.60 

-6.74 

10.00 

-0.674 

1550.00 

101.35 

-5.32 

10.00 

-0.532 

1575.00 

101,10 

-3.90 

10.00 

-0.390 

1600.00 

100.85 

-2.48 

10.00 

-0.248 

1625.00 

100.60 

-1.06 

10.00 

-0.106 

1650.00 

100.36 

0.36 

10  00 

0.036 

1675.00 

100.11 

1.78 

lv.oo 

0.178 

1700.00 

99.86 

3.19 

10.00 

0.319 

1725.00 

99.61 

4.61 

10.00 

0.461 

1750.00 

99.37 

6.03 

10.00 

0.603 

1775.00 

99.12 

7.45 

10.00 

0.745 

1800.00 

98.87 

8.87 

10.00 

0.887 

1825.00 

98.62 

8.62 

10.00 

0.862 

1850.00 

98.39 

8.39 

10.00 

0.839 

1875.00 

98.15 

8.15 

10.00 

0.815 

1900,00 

97.93 

7.93 

10.00 

0.793 

1925.00 

97.71 

7.71 

10.00 

0.771 

1975.00 

97.29 

7.29 

10.00 

0.729 

2025.00 

96.89 

6,89 

10.00 

0.689 

2075.00 

96.51 

6.51 

10.00 

0.651 

10.0  0.20000E-03 
10.0  0.20000E-03 
10.0  0.20000E-03 
10.0  0.20000E-03 
10.0  0.20000E-03 
10.0  0.20000E-03 

10.0  0.20000E-03 
10.0  0.20000E-03 
10.0  0.20000E-03 
10.0  0.20000E-03 
10.0  0.20000E-03 
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2125.00 

96.15 

6.15 

10.00 

0.615 

2225.00 

95.49 

5.49 

10.00 

0.549 

2325.00 

94.90 

4.90 

10.00 

0.490 

2425.00 

94.37 

4.37 

10.00 

0.437 

2625.00 

93.48 

3.48 

10.00 

0.348 

2825.00 

92.76 

2.76 

10.00 

0.276 

3325.00 

91.56 

1.56 

10.00 

0.156 

3825.00 

90.85 

0.85 

10.00 

0.085 

4325.00 

90.41 

0.41 

10.00 

0.041 

4825.00 

90.10 

0.10 

10.00 

0.010 

5000.00 

90.00 

0.00 

10.00 

0.000 
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APPENDIX  D:  EXAMPLE  DATA  FILES 


1.  Following  are  listings  of  "standard"  data  files  used  as  the  basis  of 
parametric  studies  in  Part  III  of  the  report.  Working  from  these  files,  the 
parametric  studies  were  performed  by  systematically  altering  the  input  file. 
Copies  of  screen  display  illustrating  the  geometry  are  shown  in  Part  III  of 
this  report. 

2.  DATACHK  models  a  levee  section  amenable  to  conventional  analysis.  It 
has  a  10-ft-thick  top  blanket  overlying  an  80-ft-thick  pervious  substratum. 

An  example  is  listed  below: 

IRREGULAR  FOUNDATION 
TEST  PROBLEM 
.2000 

2  "CONST"  .0002  110 
0  0  80  90 
1550  0  80  90 
2  "CONST"  .0002 
1800  0  80  90  90 
5000  0  80  90  90 
NO  WELL 

3.  File  DATAIRR  models  an  irregular  top  blanket  having  a  thick  clay  plug 
paralleling  the  levee.  An  example  of  such  is  listed  below: 

IRREGULAR  FOUNDATION 
TEST  PROBLEM 
0.200 

2  "CONST"  .0002  175 

750.0  60.0  140.0  158.0 

1750.0  60.0  140.0  160.0 

4  "CONST"  .0002 

1900.0  60.0  140.0  160.0  160.0 

2400.0  60.0  120.0  158.33  158.33 

2800.0  62.0  140.0  155.0  158.33 
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4900.0  70.0  140.0  158.0  158.33 
NO  WELLS 

4.  File  DATADCH  models  a  levee  section  with  a  15-ft-thick  top  blanket 
and  a  landside  ditch.  It  is  listed  below: 

TEST  FILE  "DATADCH" 

DITCH  IV  ON  3H  200  FT  FROM  LEVEE  10  FT  DEEP 
200 

2  "CONST"  .0003  100 
0  0  65  80 
2000  0  65  80 
6  "CONST"  .0003 
2200  0  65  80  80 
2400  0  65  80  80 
2430  0  65  70  80 
2440  0  65  70  80 
2470  0  65  80  80 
7000  0  65  80  80 
NO  WELLS 

5.  File  DATAWELL  is  similar  to  file  DATACHK  with  the  exception  that  a 
line  of  relief  wells  is  molded  at  the  levee  toe.  It  is  also  listed  below: 

IRREGULAR  FOUNDATION 
TEST  PROBLEM 
0.2000 

2  "CONST"  0.0002  110 
0  0  80  90 
1500  0  80  90 

2  "const"  o.ooo; 

1800  0  80  90  90 
5000  0  80  90  90 
WELL 
1820  94 
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APPENDIX  E:  HAND  CHECK,  FILE  DATACHK 
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APPENDIX  F:  PROGRAM  LISTING 


PROGRAM  LEVEENSU 

LEVEE  UNDERSEEPAGE  ANALYSIS 
FOR 

IRREGULAR  FOUNDATION  CONDITIONS 

Thoaas  F.  Holff ,  Ph.D.,  P.E. 
Michigan  State  University 
Hay  1989 

Last  revision  6/1/89 


'  tU  START  ttt 
start: 

CLEAR 

set  diaensions  — 

REN  ISTATIC 

DIN  x(40),  yl(40),  y2(40),  y3(40),  y*ater(40) 

DIN  xx(200|,  yyl(200),  yy2(200),  yy3(200),  yy«ater{200) 

DIN  d (200) ,  2(200),  kb(200),  pie2el(200),  oldpzel{200) 

DIM  cl (200) ,  c2(200),  c3(200) 

DIN  reshead(200),  iexit(200) 

‘  —  define  foreat  strings  — 

ForaStrll  =  'Iteration  ttll  i  aax  =  II. Ill  at  x  =  Hill* 

ForaStr2t  =  *  h  aax  =  II. Ill  at  x  =  Hill* 

ForaStr3t  =  'ltlll.lt  lllt.lt  till. II  till. II  till. II  III. I  III. I  t.lllltAAAA' 
ForaStrll  «  ‘Hill. II  till. II  lilt. II  till. II  III. Ill  ' 
blank40t  =  ' 

* —  set  changedate  — 
changedatel  =  '  Hay  31,  1989* 

*  —  set  defaults  and  initialize  variables  — 
dx  =  25 

aaxitersl  =  1000 
tol  =  .0005 
Mel  II  =  0 
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'  m  OPENING  SCRECH  It* 


'  —  print  banner  — 

SCREEN  0 
CIS 

COLOR  15,  2 

PRINT  *  —  LEVEENSU  — 

PRINT  *  UNDERSEEPASE  ANALYSIS 

PRINT  *  FOR 

PRINT  *  IRREGULAR  FOUNDATION  CONDITIONS  * 

PRINT 

PRINT  ‘  Thoeas  F.  Nolff 

PRINT  *  Michigan  State  University 

PRINT 

COLOR  15,  4 
PRINT 

PRINT  *  U.  S.  Arey  Engineer 
PRINT  *  Naterways  Experiment  Station 
PRINT  ’  Geotechnical  Laboratory 

PRINT 
PRINT 

PRINT  “Last  revision  *;  changedate* 

60SUB  asuflag 
G0SU8  corpsflag 


' —  set  graphics  chip  — 

COLOR  15,  2 
LOCATE  17,  1 
setchip: 

PRINT 

PRINT  'Enter  graphics  «ode* 

PRINT  'EGAC  for  Hi  Res  color' 

PRINT  'EGAN  for  Hi  Res  eonochrote' 

PRINT  "C6A  for  Ned  Res  «onochro«e‘ 

INPUT  chip* 

IF  chip*  =  'EGAC'  OR  chip*  =  'egac'  THEN 
chip*  =  ’egac* 
sandcolor  *  8 
daycolor  *  6 
leveecolor  *  4 
matercolor  *  3 
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ELSEIF  chip!  =  'EGAN*  OR  chip!  =  'egaa'  THEN 
chip!  *  *egaa’ 

Ron1(  =  CHR$ (0)  ♦  CHRI(O)  ♦  CHR$(0)  ♦  CHR$(0) 

Ron2$  =  CHR$(0)  ♦  CHRt ( 0 )  ♦  CHR$(0)  ♦  CHR$ ( 0 ) 

Rom3$  *  CHR$(16)  ♦  CHR${U)  ♦  CHR*(16)  ♦  CHRt(O) 

Ro*4$  *  CHR$ (0 )  *  CHR$ (0)  ♦  CHR${0)  ♦  CHRS(O) 

Rom5$  =  CHRi ( 0)  +  CHR$(0)  *  CHRI(O)  ♦  CHR«(0) 

sandtilel  =  Rowl$  *  Rom2$  ♦  Ro*3$  ♦  Ro*4$  +  Ro*5$ 

sandcolor  =  7 

daycolor  =  7 

leveecolor  =  7 

aatercolor  =  7 

ELSEIF  chips  ^  'CGA*  OR  chip*  =  'cga*  THEN 
chips  =  *cga* 
sandcolor  =  1 

sandtileS  =  CHRS(O)  +  CHRS(O)  ♦  CHRI(16)  ♦  CHRS{0)  +  CHRS{0) 
daycolor  =  1 
leveecolor  =  1 
ttatercolor  =  1 
ELSE 

GOTO  setchip 
END  IF 


'  III  SECOND  SCREEN  -  INPUT  III 
getinfilenaae: 

COLOR  7,  1 

CLS 

PRINT 

COLOR  10,  1 

PRINT  'ENTER  INPUT  FILE  NAME" 

COLOR  7,  1 

ON  ERROR  60T0  getinfilenaae 
INPUT  infileS 

OPEN  infilel  FOR  INPUT  AS  II 
CLOSE  1 

getoutfilenaae: 

ON  ERROR  GOTO  getoutfilenaae 
PRINT 

COLOR  10,  1 

PRINT  r ENTER  OUTPUT  FILE  NAME* 

COLOR  7,  1 
INPUT  outfileS 

IF  (outfile*  *  infilel)  THEN 
COLOR  10,  1 

PRINT  'BUST  BE  DIFFERENT  THAN  INPUT  FILE1 
60T0  getoutfilenaae 
END  IF 
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•  lit  THIRD  SCREEN  --  DEFAULTS  III 
defaults: 

COLOR  7,  1 
CIS 

COLOR  10,  1 

PRINT  'DEFAULT  SETTINGS' 

COLOR  7,  1 

PRINT  'Press  return  if  OK* 

PRINT  'Type  nee  nueber  to  change* 

PRINT 

pointl: 

ON  ERROR  GOTO  pointl 
PRINT  'Tolerance  =  ';  tol 
INPUT  value 
IF  value  >  0  THEN 
tol  =  value 
COLOR  10,  1 
GOTO  pointl 
END  IF 
COLOR  7,  1 

point2: 

ON  ERROR  GOTO  point2 

PRINT  'Haxieue  iterations  =  ':  eaxitersZ 
INPUT  value 
IF  value  >  0  THEN 

■axitersZ  =  value 
COLOR  10,  1 
60T0  point2 
END  IF 
COLOR  7,  1 

points: 

ON  ERROR  GOTO  points 

PRINT  'Naxieua  node  spacing  near  levee  =  ';  dx 
INPUT  value 
IF  value  >  0  THEN 
dx  =  value 
COLOR  10,  1 
60TO  points 
END  IF 
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•  til  III  INPUT  FROM  FILE  III  III 
readdata: 

ON  ERROR  GOTO  infileerror 
OPEN  infilet  FOR  INPUT  AS  II 
INPUT  II,  titlelt 
INPUT  II,  title2t 
INPUT  II,  kf 

INPUT  II,  nrivsecsl,  pereflagrt,  pereriv,  yriv 
FOR  jl  =  1  TO  nrivsecsl 

INPUT  II,  x(jl),  yl(jZ),  y2( jX) ,  y3(jX) 
yxater(jZ)  =  yriv 
NEXT  jX 

INPUT  II,  nlandsecsl,  penflaqlt,  peraland 
nsecsX  =  nrivsecsX  ♦  nlandsecsX 

FOR  jX  *  nrivsecsX  ♦  1  TO  nsecsX 

INPUT  II,  x(jX),  yl(jX),  y2(jX),  y3(jX),  yeater(jl) 
NEXT  jX 

INPUT  II,  nellflagt 

IF  Hellflagt  =  •well-  OR  nellflagt  =  'NELL*  THEN 
INPUT  ti,  XNell,  yvell 
wellX  *  1 
ENO  IF 
CLOSE  1 

* —  disble  error  handler  — 

ON  ERROR  GOTO  0 

'  III  FOURTH  SCREEN  -  DISPLAY  INPUT  till 
listfile: 

COLOR  7,  1 
CLS 

COLOR  10,  1 
PRINT  'Data  file: 

COLOR  7,  1 
PRINT  infilet 
PRINT 
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PRINT  titlel* 

PRINT  title2* 

PRINT 

COLOR  10,  1 

PRINT  'Foundation  Permeability* 

COLOR  7,  1 
PRINT  kf 
PRINT 

COLOR  10,  1 
PRINT  'Riverside* 

PRINT  'sections  pereflag  per*  yriv* 

COLOR  7,  1 

PRINT  nrivsecsl,  pereflagri,  pereriv,  yriv 
COLOR  10,  1 

PRINT  *x  yl  y2  y3  ' 

COLOR  7,  1 

FOR  jX  =  1  TO  nrivsecsX 
IF  CSRLIN  >  23  THEN 
COLOR  10,  1 

INPUT  'Press  return  to  continue*;  code! 

CLS 

COLOR  7,  1 
END  IF 

PRINT  x(jX),  yl(jX),  y2(jX),  y3(jX) 

NEXT  JX 
PRINT 

IF  CSRLIN  >  23  THEN 
COLOR  10,  1 

INPUT  'Press  return  to  continue';  code! 

CLS 

COLOR  7,  1 
END  IF 
COLOR  10,  1 
PRINT  'Landside' 

PRINT  'sections  pereflag  pere' 

COLOR  7,  1 

PRINT  nlandsecsX,  pereflaglt,  peraland 
COLOR  10,  1 

PRINT  '  x  yl  y2  y3 

COLOR  7,  1 
COLOR  7,  1 

FOR  jl  *  nrivsecsX  ♦  1  TO  nsecsX 
IF  CSRLIN  >  23  THEN 
COLOR  10,  1 

INPUT  'Press  return  to  continue*;  code* 

CLS 

COLOR  7,  1 
END  IF 

PRINT  x(jX),  yl(jX),  y2(jX),  y3(jX),  yvater(jl) 

NEXT  JX 
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yuater' 


IF  CSRIIN  >  23  THEN 
COLOR  10,  1 


INPUT  'Press  return  to  continue* j  code! 
CLS 

COLOR  7,  1 
END  IF 

PRINT  'WELL',  'YNELL* 

PRINT  xuell,  yuell 
END  IF 
PRINT 

COLOR  10,  1 

PRINT  'Press  Return  to  Continue* 

INPUT  code* 

*  —  FIFTH  SCREEN  -START  6RAPHICS  — 

CLS 

IF  chip*  =  'egac*  THEN 
SCREEN  9  v 

VIEN  (10,  10H630,  260),  ,  8 
ELSE IF  chip*  *  'egaa'  THEN 
SCREEN  9 

VIEN  (10,  10H630,  260),  ,  7 
ELSE 
SCREEN  2 

VIEN  (10,  10H630,  150),  ,  1 
END  IF 

' —  define  initial  uindou  — 

xuain  «  x(l)  -  100 

xuaax  =  x(nsecsl)  +  100 

yuain  =  yl(l)  -  40 

yuaax  =  y3(nsecsl)  ♦  50 

NINDON  (xuain,  yuainMxuaax,  yuaax) 

GOSUB  drausection 
GOSUB  displaytitles 


' —  pause  for  copy  — 

LOCATE  20,  1 

PRINT  'Press  return  to  continue 
INPUT  ;  code* 

‘  tit  Ml  ANALYSIS  tit  ttt 


F-7 


'  —  generate  nodes  — 
LOCATE  20,  1 
PRINT  ‘Generating  Nodes 
il  =  0 
jX  =  0 


NeaSegaent; 
iX  =  il  M 
jX  =  jZ  ♦  1 
xx(jX)  =  x(iX) 
yyl(jX)  =  yl(iX) 

yy2(jX)  =  y2(iX) 
yy3(jl)  =  y3(iX) 
yyaater(jZ)  =  yaater(iX) 
d(jX)  =  yy2(jX)  -  yyl(jX) 
z(jl)  =  yy3UX)  -  yy2(jX) 

IF  iX  <  nrivsecsX  THEN 
kb(jZ)  *  perariv 
ELSEIF  iX  =  nrivsecsX  THEN 
kb(jX)  =  perariv 
jrtoeX  =  jX 

ELSEIF  iX  =  nrivsecsX  ♦  1  THEN 
kb(jX)  *  peralanii 
jltoeX  =  jX 

ELSEIF  il  )  nrivsecsZ  +  1  THEN 
kb(jX)  *  peraland 
END  IF 


'  —  calculate  slopes  for  node  interpolation  — 

denoa  =  (x(iX  ♦  1)  -  x(iZ)) 

slopel  ■  (yl{iZ  ♦  1)  -  yl(il))  /  denoa 

slope2  =  (y2(iZ  ♦  1)  -  y2(iX) J  /  denoa 

slope3  =  (y3(iZ  ♦  1)  -  y3(iZ))  /  denoa 

slopeaater  =  (yaaterfiX  +  1)  -  ywater (iX) )  /  denoa 

NeaNode: 
jX  =  jX  ♦  1 

*  —  shorten  node  spacing  near  levee  — 

IF  xx(jZ  -  1)  <  x (nrivsecsX)  -  1000  THEN 
dxx  -  dx  t  20 

ELSEIF  xx(jZ  -  1)  <  x(nrivsecsZ)  -  400  THEN 
dxx  *  dx  t  8 

ELSEIF  xx(JZ  -  1)  <  x(nrivsecsZ)  -  300  THEN 
dxx  -  dx  t  4 

ELSEIF  xx(jZ  -  1)  <  x (nrivsecsZ)  -  100  THEN 
dxx  =  dx  t  2 

ELSEIF  xx(jZ  -  1)  >  x(nrivsecsZ  ♦  1)  ♦  1000  THEN 
dxx  *  dx  I  20 

ELSEIF  xx(jZ  -  1)  >  x (nrivsecsZ  ♦  1)  ♦  400  THEN 
dxx  =  dx  I  8 
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ELSEIF 

xx(jZ  -  1) 

> 

x(nrivsecsZ  ♦ 

1) 

♦  300 

THEN 

dxx 

=  dx  1  4 

ELSEIF 

xx(jZ  -  1) 

> 

xfnrivsecsX  + 

1) 

♦  100 

THEN 

dxx 

a  dx  1  2 

ELSE  dxx  =  dx 

END  IF 

xx(jX) 

=  XX(jZ  - 

1) 

♦  dxx 

if  ji : 

>  200  THEN 

PRINT  ‘too  aany  nodes" 

STOP 

ELSEIF  xx(jZ)  >=  x(nsecsZ)  THEN 
GOTO  StopNodes 

ELSEIF  xx(jl)  >=  x(il  ♦  1)  THEN 
JI  =  jl  -  1 
GOTO  NeuSegaent 

END  IF 

yyl(jl)  =  yyl(jZ  -  1)  ♦  slopel  I  dxx 

yy2(jZ)  =  yy2(jX  -  1)  ♦  slope2  I  dxx 

yy3(jX)  =  yy3( jX  -  1)  ♦  slope3  I  dxx 

yyaater(jZ)  =  yy*ater(jZ  -  1)  +  slopeaater  I  dxx 

d(jX)  =  yy2(jl)  -  yyl(jZ) 

z(jZ)  a  yy3(jZ)  -  yy2(jZ) 

kb(jZ)  =  kb(jZ  -  1) 

60T0  NeaNode 

StopNodes: 
nnodesZ  -  jZ 
xx(jZ)  »  x(nsecsZ) 
yyl(jX)  =  yl(nsecsZ) 

yy2(jZ)  =  y2(nsecsZ) 

yy3(jZ)  =  y3(nsecsZ) 

yyaater(jX)  =  yuater(nsecsZ) 
d(jX)  =  yy2(jZ)  -  yyl(jZ) 
z(jX)  =  yy3(jZ)  -  yy2(jZ) 
kb(jZ)  =  kb(jZ  -  1) 

LOCATE  20,  26 

PRINT  nnodesZ;  *  nodes  * 

60SUB  graphnodes 

• —  zero  blanket  peraeabilities  under  levee  — 

FOR  jZ  »  jrtoeZ  ♦  1  TO  jltoeX  -  1 
kb(jZ)  =  0 
NEXT  jZ 

• —  adjust  variable  peraeabilities  — 

IF  (peraflagrt  *  "curve")  OR  (peraflagrt  =  "CURVE")  THEN 
FOR  jZ  =  1  TO  jrtoeZ 

kb(jZ)  «  perariv  /  (EIP(-.065924  I  UO  -  z(jZ)))> 
NEXT  jZ 
END  IF 
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IF  (peraflagl*  =  ‘curve’)  OR  (peraflagl*  =  ’CURVE’)  THEN 
FOR  jX  *  jltoel  TO  nnodesX 

Itb(jX)  =  peraland  /  (EXP(-.065?24  I  (10  -  z(jZ)))) 

NEXT  jX 
END  IF 

•  lit  PRINT  NODE  INFO  TO  OUTPUT  FILE  ttt 

OPEN  outfile*  FOR  OUTPUT  ftS  12 

PRINT  12,  ’PROGRAM  LEVEEMSU  changedatet;  ’  edition’ 

PRINT  «2,  ’INPUT  FILE:  ’;  infile* 

PRINT  12,  ’OUTPUT  FILE:  ’;  outfile* 

PRINT  12, 

PRINT  12,  titlel* 

PRINT  12,  title2* 

PRINT  12, 

PRINT  12,  ’KF  =  ’;  kf 

PRINT  12,  ’PERMFLAGR  =  ’;  peraflagr*;  ’  PERHRIV  =  *}  perariv 

PRINT  12,  ’PERMFLA6L  =  ’;  peraflagl*;  *  PERMLAND  =  *;  peraland 

IF  wel 1Z  =  I  THEN 
PRINT  12, 

PRINT  12,  ’NELL  LINE  AT  X  *  ’;  xaell 

PRINT  12,  ’AVGERAGE  PIEZ  EL  =  *;  yuell 
END  IF 
PRINT  12, 

PRINT  *2,  ’  xx  yyl  yy2  yy3  yyaater  d  z  kb’ 

PRINT  12, 

FOR  jX  =  1  TO  nnodesX 

PRINT  12,  USING  ForaStr3*;  xx(jl);  yyl(jZ);  yy2(jZ)j  yy3(jl);  yyuater(jZ);  d(jl);  z(jl);  kb(jZ) 
NEXT  jZ 

FIND  NELL  NODE  IF  PRESENT  — 

IF  wel IX  =  1  THEN 

FOR  jX  =  jltoel  TO  nnodesX 

IF  (xwel  1  >  xx(iX  -  D)  AND  (xwell  <  xx(jZ  ♦  D)  THEN 
j Mel  IX  =  jl 
END  IF 
NEXT  jX 
END  IF 

INITIALIZE  PIEZOMETRIC  ELEVATIONS 
LOCATE  20,  1 

PRINT  ’Initializing  ’ 

xl  =  SQR((kf  /  kbljrtoeX  -  1))  I  z(jrtoel)  I  d(jrtoeXJ) 
cr  =  1  /  xl 

x3  =  SQR( (kf  /  kb( jltoeX  ♦  1))  t  z(jltoeZ)  I  d(jltoeX)) 
cl  =  1  /  x3 

x2  a  xx(jltoeZ)  -  xx(jrtoeZ) 
s  *  xl  ♦  x2 

h  *  yyMater(jrtoeX)  -  yyMater( jltoeX) 

a  =  h  /  (s  ♦  x3) 

hO  =  (hi  x3)  /  (s  ♦  x3) 

hOr  *  (h  I  xl)  /  (s  ♦  x3) 


F- 10 


piezel(l)  :  ywater(l) 
oldpzel (1)  =  ynater(l) 

FOR  jX  8  2  TO  nnodesl  -  1 
IF  NellX  8  1  AND  jX  -  jnellX  THEN 
piezel(jX)  -  ywell 
ELSEIF  xx(jX)  <=  xx(jrtoeX)  THEN 
xr  *  xx(jrtoeX)  -  xx(jX) 

piezel(jX)  =  yyxater( jrtoeZ)  -  hOr  I  EXP(-1  t  cr  t  xr) 

ELSEIF  xx(jX)  >=  xx(jltoeX)  THEN 
xl  =  xx(jX)  -  xx(jltoeX) 

piezel(jX)  8  yynater(jltoel)  ♦  hO  I  EXP(-1  I  cl  I  xl) 

ELSE 

piezel(jX)  8  yyuater(jrtoel)  -  hOr  -  «  I  (xx(jl)  -  xx(jrtoeX)) 
kb(jX)  8  0 
END  IF 

oldpzel(jX)  =  piezel(jX) 

NEXT  jX 

piezel(nnodesX)  =  ywater (nsecsX ) 
oldpzel (nnodesl)  8  yxater(nsecsX) 

'  —  calculate  node  constants  — 

FOR  jX  8  2  TO  (nnodesl  -  1) 

d(jX)  8  kf  I  (d(jl)  ♦  d(jX  -  1))  I  .5  /  (xx(jX)  -  xx(jX  -  1)1 
c2(jl)  8  kf  I  (d(jl  ♦  1)  +  d( jX) )  I  .5  /  ( xx ( jX  ♦  1)  -  xx(jl)) 

IF  jX  8  2  THEN 

c3(jX]  *  kb(jX)  t  ( ( ( xx (3)  -  xx(2))  I  .5)  ♦  (xx(2)  -  xx(l)I)  /  z(2) 

ELSEIF  jl  8  nnodesl  -  1  THEN 

c3(jX)  8  kb(jX)  I  ( ( ( xx( jX  ♦  l)  -  xx(jX))  ♦  (xx(jl)  -  xx(jX  -  1))  I  .5))  /  z(jl) 
ELSEIF  jl  8  jrtoeX  THEN 

c3(jl)  8  kb(jl)  t  (xx(jl)  -  xx(jX  -  1))  I  .5  /  z(jX) 

ELSEIF  jl  8  jltoeX  THEN 

c3(jX)  8  kb(jX)  t  ( xx ( jX  ♦  1)  -  xx(jX))  I  .5  /  z(jl) 

ELSE 

c3(jl)  *  kb(jX)  »  (**(:!  *  »)  -  xx(jX  -  1))  t  .5  /  z(jX! 

END  IF 
NEXT  jl 


'tit  ttt  SOLVE  ttt  ttt 
iter!  8  0 
LOCATE  20,  1 
PRINT  'Solving 

solve: 

iexiteax  *  0 
taxres  *  0 

‘ —  iterate  forward  - 
iter!  8  iter*  ♦  1 
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FOR  jl  *  2  TO  nnodesl  -  1 

piezel(jl)  =  (piezel ( jl  -  1)}  I  d(jl)  +  (piezelljl  ♦  1))  I  c2(jl)  ♦  (yynater(jZ))  I  c3(jl) 
piezel ( jZ)  =  piezel(jZ)  /  (cl(jZ)  ♦  c2(jl)  ♦  c3(jl|) 
piezel ( jZ)  =  piezel (jZ)  ♦  .2  I  (piezel ( jZ}  -  oldpzeK jl) ) :  REM  overshoot 
IF  Neill  --  1  AND  jl  =  jwellX  THEN 
piezel ( jX)  =  yvell 
END  IF 
NEXT  jl 

' —  iterate  backaard  — 

iterl  =  iterX  ♦  1 

FOR  jX  *  nnodesl  -  1  TO  2  STEP  -1 

piezelljl)  =  (piezelljl  -  1))  I  cl( jX)  +  (piezelljl  ♦  1))  I  c2(jl)  +  (yywaterl jZ) )  I  c3( jZ) 
piezelljl)  =  piezelljl)  /  (cl(jl)  ♦  c2(j!)  ♦  c3(jl)( 

IF  Nelli  =  1  AND  jl  =  jwellZ  THEN 
piezel(jZ)  =  y«ell 
END  IF 

iexit(jZ)  =  (piezelljl)  -  yyNater(jl))  /  (z(jl)) 

IF  iexitijl)  >  iexitiax  THEN 
iexitaax  =  iexit(jl) 
xiexitaax  =  xx(jl) 

END  IF 

res  =  ABS(piezelljl)  -  oldpzel ( jZ) ) 

IF  res  >  aaxres  THEN 
•axres  =  res 
xaaxres  =  xx(jZ) 

END  IF 
NEXT  jl 

set  old  piez  els  to  nea  ones 
FOR  jZ  =  2  TO  nnodesl  -  1 
oldpzel! jl)  s  piezelljl) 

NEXT  jl 

'  —  print  iteration  and  aax  gradient  — 

IF  iterl  /  10  =  FlXliterl  /  10)  THEN 
LOCATE  21,  1 

PRINT  USIN6  ForaStrlt j  iterl,  iexitaax,  xiexitaax 
END  IF 

IF  iterl  >=  aaxitersl  THEN 
GOTO  stopiter 
ELSE IF  aaxres  >  tol  THEN 
GOTO  solve 
END  IF 
stopiter: 
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'  —  calculate  residual  heads  -  - 

■axreshead  =  0 

FOR  jX  =  1  TO  nnodesl 

reshead(jl)  =  pie.eujl)  -  yynater(jX) 
IF  reshead(jX)  >  aaxreshead  THEN 
■axreshead  =  reshead(jl) 
xaaxreshead  -  xx(jX) 

FND  IF 
NE*f  jX 


' —  calculate  well  flou  — 

IF  aellZ  =  1  THEN 
jX  *  jnellX 

Qin  =  kf  t  ((piezeKjX  -  1)  -  piezel(jX))  /  (xx(jX)  -  xx(jl  -  1)1)  t  (d(jl)  ♦  d(jX  -  1))  /  2 
Qout  =  kf  t  ( ( piezel ( jX)  -  piezeKjX  ♦  1))  /  (xx(jX  ♦  1)  -  xx(jl)))  I  (d(jl  +  1)  ♦  d(jX)!  /  2 
Qup  =  kb(jX)  t  ((piezel(jX)  -  yywater ( jl) )  /  z(jX))  I  (xx(jl  +  1)  -  xx(jl  -  11)  /  2 
Quell  -  Qin  -  Qout  -  Qup 
Qugal  -  Quell  I  7.48 
END  IF 

lit  WRITE  FINAL  HEADS  AND  GRADIENTS  TO  FILE  till 
PRINT  12, 

PRINT  12,  *  xx  piezel  reshead  z  i* 

PRINT  12, 

FOR  jX  =  1  TO  nnodesl 

PRINT  12,  USING  For«Str44;  xx(jl);  piezelljX);  resheadfjX);  zfjX);  iexit(jl) 

NEXT  jX 

IF  uellX  =  1  THEN 
PRINT  12, 

PRINT  12,  USING  'Nell  flon  =  #*1.111  ftA3/ain/ft’;  Ouell 
PRINT  #2,  USING  ‘Nell  flon  =  ###.###  gpa/ft’;  Qngal 

END  IF 
CLOSE  2 

'  III  DISPLAY  FINAL  RESULTS  III 
LOCATE  20,  1 

PRINT  ’SOLUTION  COMPLETE  ’ 

LOCATE  20,  50 

PRINT  ’OUTPUT  FILE:  outfile# 

LOCATE  21,  1 

PRINT  USING  ForaStrl*',  iterX,  iexitaax,  xiexitaax 
LOCATE  22,  1 

PRINT  USING  ForaStr2*;  aaxreshead,  xaaxreshead 
IF  nellX  =  1  THEN 
LOCATE  21,  48 

PRINT  USING  ’Nell  flon  =  III. Ill  gpa/ft’j  Qugal 
END  IF 

60SUB  graphpiez 
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*  —  check  for  nett  window  — 
checkwindow: 

LOCATE  23,  1 

PRINT  blanket;  blanket; 

LOCATE  24,  1 

PRINT  blank40*;  blank40*; 

LOCATE  23,  1 

INPUT  ;  'Change  graphics  window  (Y  or  N)  code* 

IF  code*  =  'V  OR  code*  =  y  THEN 
LOCATE  23,  1 

PRINT  'Current  xwain,  xwaax,  ywain,  ywaax:*; 

LOCATE  23,  40 

PRINT  USING  '  «!«««  «««««  *«•»«  HIM';  xwain,  xwaax,  ywain,  ywaax; 
LOCATE  24,  1 

PRINT  ‘Enter  new  values  with  coaaas:  '; 

LOCATE  24,  40 

INPUT  ;  xwain,  xwaax,  ywain,  ywaax 
CLS 

MINDOH  (xwain,  ywain)-(xwaax,  ywaax) 

60SU8  drawsection 
60SUB  graphnodes 
SOSUB  graphpiez 
GOSUB  displaytitles 
GOTO  checkwindow 
END  IF 

‘ —  print  output  file  to  printer  if  desired  — 
printer: 

LOCATE  23,  1 

INPUT  ;  'Duap  output  to  printer  (Y  or  N)  ';  code* 

IF  code*  =  'Y*  OR  code*  =  'y'  THEN 
GOTO  printer2 

ELSE1F  code*  =  'N'  OR  code*  =  'n'  THEN 
GOTO  recycle 

ELSE  60T0  printer 
END  IF 

printer2: 

LPRINT  "PROGRAM  LEVEEMSU  ';  changedate*;  •  edition' 

LPRINT  DATE* 

LPRINT  TINE* 

LPRINT  'INPUT  FILE:  ';  infile* 

LPRINT  'OUTPUT  FILE:  ';  outfile* 

LPRINT 

LPRINT  titlel* 

LPRINT  title2* 

LPRINT 

LPRINT  *KF  *  kf 

LPRINT  'PERHFLAGR  *  ';  peraflagr*;  ■  PERHRIV  «  ';  perariv 
LPRINT  'PERHFLAGL  *  ';  peraflaglt;  ■  PERNLAND  *  ';  peraland 
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IF  Neill  *  1  THEN 
LPRINT 

LPRINT  ’NELL  LINE  AT  I  *  'j  xeell 
LPRINT  'AVGERA6E  PIEZ  EL  *  *5  yeell 
END  IF 
LPRINT 

LPRINT  *  xx  yyl  yy2  yy3  yyeater  d  z  kb* 

LPRINT 

FOR  jl  =  1  TO  nnodesl 

LPRINT  USING  ForaStr3$;  xx(jl):  yyl(jl);  yy2( jl) j  yy3( jl) j  yywater( jX) j  d(jl);  z(jl);  kb(jl) 

NEXT  jl 

LPRINT 

LPRINT  '  xx  piezel  reshead  z  i' 

LPRINT 

FOR  jX  *  1  TO  nnodesl 

LPRINT  USING  ForaStrll;  xx(jl);  piezel(jl);  reshead(jl);  z(jl);  iexit(jl) 

NEXT  jl 

IF  Neill  ;  1  THEN 
LPRINT 

LPRINT  USIN6  'Nell  floN  -  III. til  ftA3/ain/ffj  Quell 
LPRINT  USING  'Nell  TIon  =  (««.«««  gpa/ft'j  Segal 
END  IF 

'  —  check  for  nee  problea  — 
recycle: 

LOCATE  24,  1 

INPUT  ;  'Run  again  (Y  or  N)  '}  code! 

IF  code!  =  T  OR  code!  =  'y*  THEN 
GOTO  start 

ELSEIF  code!  =  'N*  OR  code!  =  V  THEN 
GOTO  terainate 
ELSE  GOTO  recycle 
END  IF 

terainate: 

END 


- SUBROUTINES - 

inf ileerror : 

PRINT  'error  in  input  file' 


errorhandler: 

PRINT  'fatal  error 


*  —  HSU  Flag— 

■suf lag: 

COLOR  IS,  2 
FOR  il  =  1  TO  7 
LOCATE  il,  50 
PRINT  * 

NEXT  il 
xf  =  CHRt(2I9) 

LOCATE  2,  55 

PRINT  xt;  xt;  xt;  xS;  xt 
LOCATE  3,  55 
PRINT  xt 
LOCATE  4,  55 

PRINT  xt;  xt;  xt;  xt;  xt; 

LOCATE  5,  59 
PRINT  xt; 

LOCATE  6,  55 

PRINT  xt;  xt;  xt;  xt;  xt; 

RETURN 

'  -  CORPS  FLAG  — 
corps flag: 
st  -  CHRt(32) 

COLOR  15,  4 
LOCATE  9,  50 

PRINT  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st; 
LOCATE  10,  50 

PRINT  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st; 
LOCATE  11,  50 

PRINT  st;  st;  st;  xt;  xt;  st;  xt;  xt;  xt;  st;  xt;  xt;  st;  st;  st; 
LOCATE  12,  50 

PRINT  st;  st;  st;  st;  xt;  xt;  xt;  xt;  xt;  xt;  xt;  st;  st;  st;  st; 
LOCATE  13,  50 

PRINT  st;  st;  st;  xt;  xt;  xt;  xt;  st;  xt;  xt;  xt;  xt;  st;  st;  st; 
LOCATE  14,  50 

PRINT  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st; 
LOCATE  15,  50 

PRINT  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st;  st; 
RETURN 

‘  lit  Plot  Levee  Section  tit 
(transection: 

•  —  outline  pervious  substratu*  — 

LINE  (x(l),  yldlHxd),  y2(l)|,  sandcolor 

LINE  (x(nsecsl),  yt(nsecsZ))-(x(nsecsl),  y2(nsecsZ}},  sandcolor 

FOR  jZ  5  1  TO  nsecsZ  -  1 
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LINE  (x(JX),  yl(jX)Hx(jI  *  1),  yllil  ♦  U),  sandcolor 
LINE  (x(jl),  y2(jX)Hx(jX  ♦  1),  y2(jl  ♦  1)),  sandcolor 
NEXT  jX 

' —  paint  pervious  substratui  — 
ya  *  (y2(nrivsecsX  ♦  1)  ♦  yeain)  /  2 
yb  *  (y2(nrivsecsX  ♦  1)  ♦  yllnrivsecsX  ♦  1) )  /  2 
IF  ya  >  yb  THEN 
y s  ya 
ELSE 
y s  yb 
END  IF 

IF  chip*  =  'eqac*  THEN 

PAINT  (xlnrivsecsX  ♦  11,  y),  sandcolor 
ELSEIF  chip*  *  *eqaa’  THEN 

PAINT  (xlnrivsecsX  ♦  1),  y),  sandtile*,  7 
ELSE 

PAINT  (xlnrivsecsX  ♦  1),  y),  sandtile* 

END  IF 

‘  —  outline  iapervious  topstratua  — 

LINE  (x(l),  y2(l)Hx(l),  y3(l)l,  daycolor 

LINE  (x(nsecsX),  y2(nsecsX)Hx(nsecsX),  y3(nsecsX)),  claycolor 

FOR  jX  *  1  TO  nsecsX  -  1 

LINE  (x(jX),  y2(jX))~(x(jX  ♦  D.  y2(jX  ♦  1)1,  claycolor 
LINE  (x(jX),  y3(jl))'(x(jX  ♦  1),  y3(jX  ♦  1)1,  claycolor 
NEXT  jX 

*  —  paint  topstratua  — 

IF  chip*  =  'egac*  THEN 

FOR  jX  =  2  TO  nsecsX  -  1 
y  =  (y2(jX)  +  y3(jX))  /  2 
PAINT  (x(jX),  y),  claycolor 
NEXT  jl 
END  IF 

*  —  levee — 
xrtoe  =  x(nrivsecsl) 
xltoe  =  xlnrivsecsX  +  1) 
xlevaid  =  (xrtoe  ♦  xltoe)  /  2 
yrtoe  *  y3(nrivsecsX) 

yltoe  =  y3(nrivsecsX  ♦  1) 
ylevaid  =  (y3(nrivsecsl)  ♦  y3(nrivsecsX  ♦  11)  12 
LINE  (xrtoe,  yrtoeHxlevaid,  yriv),  leveecolor 
LINE  (xlevaid,  yriv)-(xltoe,  yltoe),  leveecolor 
LINE  (xrtoe,  yrtoe)*(xltoe,  yltoe),  leveecolor 
PAINT  (xlevaid,  (ylevaid  ♦  yriv)  /  2),  leveecolor 

*  —  eater  — 

LINE  (x(l),  yrivMxrtoe,  yriv),  aatercolor 
FOR  iX  *  nrivsecsX  ♦  1  TO  nsecsX  -  1 

LINE  { x ( iX) ,  y*ater(iI))-(x(iX  ♦  1),  ywaterlil  ♦  U),  aatercolor 
NEXT  iX 
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LOCATE  l,  39 
PRINT  ywaax ; 
LOCATE  19,  39 
PRINT  yminj 
LOCATE  12,  1 
PRINT  xaain; 

LOCATE  12,  74 
PRINT  xaaax; 

RETURN 


*  —  display  titles  — 
displaytitles: 

LOCATE  1,  1 

PRINT  'INPUT  FILE:  infilet 

LOCATE  2,  1 

PRINT  titlelt 

LOCATE  3,  1 

PRINT  title2t 

LOCATE  1,  65 

PRINT  DATES 

LOCATE  2,  65 

PRINT  TINES 

RETURN 

* —  graph  nodes  — 
graphnodes: 

radius  =  .005  t  (xaaax  -  xuain) 

aspect  =  .7 

FOR  jl  =  1  TO  nnodesZ 

CIRCLE  (xx(jZ),  yy2(jZ)),  radius,  ,  ,  ,  aspect 
NEXT  jZ 
RETURN 

* —  graph  piezoaetric  line  — 
graphpiez: 

FOR  jZ  =  1  TO  nnodesZ  -  1 

LINE  (xx(jZ),  piezell jZ) )-(*x( jZ  ♦  l),  piezelljZ  ♦  II) 
NEXT  jZ 
RETURN 

END 
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APPENDIX  G :  NOTATION 


d  Thickness  of  pervious  substratum 
h0  Residual  head  at  levee  toe 

hx  Residual  head  at  distance  x  from  landside  toe  levee 
H  Net  head  on  levee 
kb  Permeability  of  top  blanket 

kbl  Permeability  of  landside  top  blanket 

kbr  Permeability  of  riverside  top  blanket 

kf  Permeability  of  pervious  substratum 

Li  Distance  from  open  seepage  entrance  to  riverside  levee  toe 
L2  Distance  from  riverside  levee  toe  to  landside  levee  toe 
L3  Distance  from  landside  levee  toe  to  open  exit 
Q  Flow 

x3  Distance  from  landside  levee  toe  to  effective  seepage  exit 
z  Thickness  of  top  blanket 
MSL  Mean  Sea  Level 
MGL  Mean  Gulf  Level 
NGVD  National  Geodetic  Vertical  Datum 
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